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1. Introduction 

Parton shower algorithms with hadronization models provide a way of generating simu- 
lated events according to approximations based on QCD. Since complete final states are 
generated, one can generate completely exclusive cross sections in this way. By summing 
over variables that one chooses not to examine, one can also make predictions for inclusive 
observables. Of special interest are predictions that, in a perturbative expansion, involve 
large logarithms of ratios of different momentum scales in the physical problem. 

An important example is the distribution of the transverse momentum p± of Z-bosons 
produced in hadron-hadron collisions at some fixed rapidity Y, da/(dp± dY). When p\ <C 
M|, the perturbative expansion of this cross section contains two powers of the large 
logarithm log(p^/M^) per power of a s . The large logarithms spoil the usefulness of fixed 
order perturbation theory for this observable and for observables that contain similar large 
logarithms. A parton shower algorithm sums contributions to the desired cross section that 
contain arbitrarily high powers of a s ; thus it sums the accompanying logarithms. For this 
reason, one can hope that a parton shower calculation will do better than a fixed order 
perturbative calculation in the region p\ <C M§ and analogous regions for other processes. 

Indeed, the basic approximation in a parton shower is that one parton splits into two 
daughter partons with a probability that matches the singularities of the QCD matrix 
element when the two daughter partons are collinear or one of them is soft. It is just these 
soft /collinear configurations that give rise to the large logarithms. Thus one can hope that 
the cross section generated by a parton shower will be a good approximation to the true 
QCD result. 

In many cases, including the Z-boson transverse momentum distribution, there are 
predictions based on the full field theory. That is, we know that the cross section "ex- 
ponentiates" in a sense that one can state precisely and we know some of the leading 
coefficients that appear in the exponent of the formula that expresses the QCD result. 
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When, for a particular process, one knows the summation of large logarithms in full 
QCD, then it is of significant interest to investigate whether a given shower algorithm 
produces matching results. To do this, one needs to derive the corresponding summation 
in the shower model, deriving the appropriate evolution equation for the observable in 
question from the general evolution equation for the shower algorithm. 

We have argued above that, because parton shower algorithms are generally based 
on parton splitting probabilities that have the proper soft /collinear singularities, these 
algorithms may provide good approximations to the full QCD result in particular cases 
involving summing large logarithms. However, a parton shower algorithm contains several 
ingredients beyond the parton splitting probabilities. Among theses are the momentum 
mapping, the color and spin treatment, and the choice of evolution variable. Depending on 
the choice of these ingredients, one may obtain agreement with full QCD for a particular 
observable or one may fail to obtain agreement. Certainly, it is widely understood that 
a virtuality ordered parton shower without a proper inclusion of the effects of quantum 
interference can get results for some observables that do not match QCD. In contrast, 
there is an extensive body of literature [1, 2, 3, 4, 5, 6, 7, 8, 9] that suggests that parton 
showers based on ordering in parton emission angles does better for many observables. In 
this paper, we will indeed see that, for the Z-boson transverse momentum distribution, the 
choice of ingredients matters. 

Because the choice of ingredients matters, we think it important to validate shower 
evolution schemes against known results for summations of large logarithms. We believe 
that such a validation program could be useful for understanding the range of validity 
of current parton shower event generators and could provide important guidelines for the 
current and future development of such programs. 

In ref. [10], in response to ref. [11], we investigated the parton energy distribution in 
electron-positron annihilation as predicted by the virtuality ordered, color-dipole based par- 
ton shower algorithm of refs. [12, 13, 14] and by fox-ordered variants of this. We concluded 
that the predictions of these parton shower algorithms are consistent with the field theory 
result represented by the well know DGLAP evolution equation [15]. This conclusion was 
confirmed by ref. [16], which included numerical studies. The parton energy distribution is 
an example of the general program of summing large logarithms, namely the logarithm of 
the resolution scale for finding the partons (~ jets) divided by the electron-positron energy. 
However, this is a rather simple example in that there is only one logarithm per power of 
a s . Observables for which there are two logarithms per power of a s present a much more 
subtle problem. 

In this paper we investigate the transverse momentum distribution of a Z-boson pro- 
duced in hadron-hadron collisions, the Drell-Yan process. Here, there are two logarithms 
per power of a s . We use a slightly modified version 1 of the shower evolution of ref. [12]. 
This evolution equation describes the evolution of the partonic states in a fully exclusive 
way. We manipulate the shower equation to produce an evolution equation for the trans- 

*We will find that we need to change the momentum mapping for initial state radiation from that of 
ref. [12]. In addition, one of the choices given in ref. [14] for a certain function A' lk that was left unspecified 
in ref. [12] gives a satisfactory result, but another of the choices does not work. 
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verse momentum distribution of the Z-boson and we compare the result to the well known 
field theory prediction that is given in ref. [17]. We find that the virtuality ordered shower 
works well in reproducing the full QCD result. 

The shower evolution equation of ref. [12] includes quantum interference among colors 
and spins. We note that this shower evolution equation is not directly practical for gener- 
ating events. A simple approximation to this evolution equation is to average over spins 
and take the leading color approximation, which yields an evolution equation [13] that 
can be written as a Markov process and is thus directly practical for generating events. 
We treat the full evolution equation, but we will see that the same results with some 
small adjustments hold for the Z-boson transverse momentum distribution produced by 
the spin-averaged, leading-color shower evolution. 

The investigation that we carry out in this paper, and that we believe would be useful 
for other observables and other shower algorithms, is analytical. That is, we want to 
test whether the transverse momentum distribution produced by the shower algorithm 
has the proper exponential form and, if so, whether the coefficients in the exponent are 
correct. A followup study, not addressed in this paper, would be numerical: how well 
does an actual implementation of the parton shower algorithm produce results that match 
numerical results given by the QCD formula. In this case, the summed QCD results, 
including nonperturbative parameters that are fit to experiment, could be obtained from 
the Resbos code [18]. We believe that one should start with an analytical study rather 
than a numerical study for the following reason. The parton shower, for the case of the 
production a Z-boson with small p±, represents the physics on "hard" scales from Mz to 
a few GeV. This is perturbative physics that is adjustable only to a limited extent. In 
particular, one can modify the parton splitting probabilities in a fashion that does not 
change them in the soft and collinear limits. 2 A full parton shower event generator also 
contains elements, such as a hadronization model and a model for the underlying event, 
that represent soft scale physics. One can adjust the models for the soft physics by tuning 
various parameters. For very large A/|/Aq CD and M^/pj_, only the part of the shower 
algorithm that cannot be tuned should matter. However, for realistic values of these 
parameters, the tunable parts of the parton shower can make a numerical difference. One 
would not want to tune the parameters in order to repair a numerical disagreement that 
was actually caused by a mismatch between the parton shower and full QCD with respect 
to the hard scale physics that is fixed. The way to check the match of hard scale physics 
is to compare analytically. Once this is checked, the numerical comparison is appropriate 
and needed. 

We organize this paper as follows. In sections 2,3, and 4, we review briefly the needed 
features of the shower evolution of ref. [12] and set up the notation for our analysis. Then 
in section 5 we outline the derivation to come and state the nature of the approximations 
that we will need. The derivation is given in sections 6, 7, 8, 9, and 10. The solution of 
the evolution equations is given in section 11 and the result is compared to the full QCD 
result in section 12. We discuss what one would get with other sorts of shower evolution in 

2 The Vincia parton shower algorithm [19] takes advantage of this flexibility. 
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section 13 and summarize the results in section 14. In appendix A, we summarize results 
from refs. [12, 13, 14] that are used in this paper. In appendix B we prove a certain property 
of integrals of the Jo Bessel function. 

2. Shower states and shower evolution 

We analyze the transverse momentum distribution of a Z-boson as generated by the parton 
shower evolution equations of ref. [12]. The organizing principle of the shower is that one 
starts at the hard interaction q + q^Z and moves to softer interactions, always factoring 
the softer interaction from previous harder interactions. 

In the notation of ref. [12], states in the sense of statistical mechanics are repre- 
sented by ket vectors |p), while possible measurements are represented by bra vectors 
(f\. Thus (F p) is the cross section that one obtains a particular result F from a mea- 
surement on an ensemble of systems represented by p). We use basis states labelled by 
lists {p, f, s' , c' , s, c} m of parton quantum numbers for two initial state partons and m final 
state partons. As the shower progresses, m increases. The momenta of the final state 
partons are {pi, . . . ,p m }- The first final state parton is the Z-boson, with momentum 
Pi = Pz- The momenta of the initial state partons are specified by giving their momentum 
fractions, r? a and rfo. Then their momenta are 

Pa = WA , {21) 
Pb = VbPB , 

where pa and pe are the momenta of the incoming hadrons, treated as massless. The 
flavors g, u, u,d,d, . . . of the final state partons are {f\, . . . , f m }, while the flavors of the 
initial state partons are denoted by a and b. The spins are specified by 

{s , s} m = {s a , Sb, s l> • • • ) s mi s a' s b' s l> ■ • ■ i s ml ' (2-2) 

with two spin labels for each parton because we use the quantum density matrix in spin 
space in order to represent possible interference among spin states. The colors are simi- 
larly specified by {c', c} m . Quantum spin and color states are represented by bra and ket 
vectors with angle brackets, as in the quantum spin inner product ({s'} m |{s} m V The lists 
{p, /, s' , c', s, c} m include all of these parton quantum numbers. The unit operator in the 
space of statistical states can be written as 

1 = ^^7/ [ d {P'f' S '' C '> S ' C }rn] \{P,f,s',C,S,c} m )({p,f,s',C,S,c} m \ , (2.3) 

where [d{p, f, s', c', s, c} m ] indicates integrating and summing over all of the quantum 
numbers. A much more detailed specification of our notation is provided in ref. [12]. In 
ref. [12], the quarks can have non-zero masses, but in this paper we take all of the partons 
except for the Z-boson to be massless. 

Let |p(t)) = U(t, 0) p(0)) be the statistical state at shower time t. Here t = at the 
beginning of the shower, which starts with the hard process q + q — > Z, and t = £f gives 
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(2.4) 



the statistical state at the end of the shower, before hadronization. (We do not discuss 
a hadronization model in this paper.) The total cross section for producing a Z-boson is 
expressed using the vector (l|, which represents the totally inclusive measurement, 

<*Z= = I [ d {P>/) S V> s > c }m] 

x /, s', c, s, c} m ) ({p, /, s', c, s, c} m \p(t)) 

= I [ d {Pifi s ', c ', s > c }m] ({s'} m \{s} m ) ({c'} m \{c} r 

x {{pj,s',c,s,c} m \p(t)) . 

Here we use the definition 

(l|{p> /, s', c', s, c} m ) = ({s'} m \{s} m ) ({c'} m \{c} m ) . (2.5) 

The total Z-production cross section is independent of t: the shower evolution maintains 
(l|W(t,0) = (1|. 

We are interested in the differential cross section da/(dp±dY) as obtained in the 
shower at the shower final time if, 

^-(PX,^)). (2.6) 

Here, the measurement function (p±, Y\ measures the cross section for the Z-boson to have 
transverse momentum p± and rapidity Y. 3 The definition is 

{ P± ,Y\p(t))= [ [d{p,f,s',c',s,c} m ] 

X (P±,Y\{P, f, S, C, S, C} m ) ({p, f, S, C, S, C} m \p) 
= / [ d fe/' s ') c '> s ' c }m]{{s'} m |{s}m) ({ c '}m\{c}m) (2.7) 

xS( Pz +- P± )6(ho g ( P -^)-Y) 
\2 \pz-PAj J 

X {{P,f,s',C,S,c}m\p{t)) . 

The starting point of shower evolution is the Born cross section, 

da 



dp± dY 



= {p±,Y\p(0)) , (2.8) 

Born 



which is proportional to 5{p±). (See eq. (4.17) in section 4.5 for details.) 
The shower evolution is specified in ref. [12] in the form 

±\p(t)) = [Hi(t)-V(t)]\p(t)) . (2.9) 



3 In general, we denote vectors in two transverse dimensions by boldface letters like p± . The transverse 
part of a four- vector p is p±. Then p\ — —p\ < 0; however, to avoid confusion, we avoid writing p\. 
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Here TCi(t) is the splitting operator, which takes a basis state with m final state partons 
and changes it to a state with m + 1 final state partons. Next, V(t) is a "virtual splitting" 
operator that leaves number of partons and their momenta, flavors, and spins unchanged. 4 
The operator V(t) is determined from 7i\{t) in such a way that 

(l|V(t) = (l|Wi(t) • (2.10) 

With this condition, 

|(l|p(t))=0 , (2.11) 

so that the total cross section to produce a Z-boson remains the Born cross section, even 
though the Z-boson momentum changes as a result of recoils against parton splittings in 
the shower. In particular, the Z-boson transverse momentum distribution broadens as the 
shower develops from t = to t = tf . We need to follow the shower evolution to find how 
the transverse momentum distribution broadens. 



3. Initial state splitting kinematics 

We will first need some kinematics for the description of an initial state splitting. For 
notational convenience, we suppose that it is parton "a" that splits. 




Figure 1: Illustration of the kinematics for an initial state splitting. 

The initial state parton with momentum p & splits, in backward evolution, to a new 
initial state parton with momentum p & and a final state parton with momentum p m +i, as 
illustrated in figure 1. The other initial state parton has momentum p^ before the splitting 
and momentum p\, equal to p^, after the splitting. In this paper, 5 we describe the splitting 

4 In general, the operator V(t) is a non-trivial operator on the partonic color space. In the leading color 
approximation, valid for N c — > oo, it is diagonal in the color space, as described in ref. [13]. The derivation 
in this paper is given for the exact treatment of color, but works also in the leading color approximation. 

In ref. [12], the splitting operator is expressed in terms of momenta rather than splitting variables y, z, <j> 
and we did not specify a choice of splitting variables. In ref. [13], we did make a choice, but that choice is 
not the same as the choice used here. 
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using splitting variables (y,z,(f>) or, alternatively, (y,k±), denned by 

1 — z . y , 

p m +i = (l + y)p a , + z — — p h + k± , 

z l + y 

i + y (3.1) 

Pa = Pa, , V ^ 

Z 

Pb = Pb ■ 

Here /cj_ is the part of p m +i that is orthogonal to both p a and pt>- 
Virtuality. The variable y is 

2p m+r p a 

y = -s • (3.2) 



We will use a shower time i based on virtuality. 



t = \og(— , 3.3 

where M is the Z-boson mass. We will typically use t as our virtuality variable instead of 
y, so that y is 

M 2 



V=^^e- t . (3.4) 



(See also eq. (3.23) below. 



Momentum fraction. The fraction of the momentum p a in the direction of p a that is 
carried away by the emitted final state parton m + 1 is 

(3.5) 

PaPb 

The variable z must be in the range < z < 1. The momentum fraction r/ a of parton "a" 
has a new value after the splitting. From eqs. (2.1) and (3.1), we have, using 2p a -pb = V&VbS 
with s = 2pa-Pb, 



f?a = T? a = - 

z z 
fib = lib ■ 



M 2 _, 

?? a + e 

rjbS 



(3.6) 



These are the exact relations. In our applications, we will generally neglect M 2 e~ t com- 
pared to rj^rjbS. That is, we consider the virtuality M 2 e~ t of a splitting to be small com- 
pared to the momentum scale of the hard process, M 2 , which, in turn, is always smaller 
than i] a i] h s. 



"For a final state splitting of parton I, the shower time is t = log(M 2 /(2p m +i -pi)). 
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Transverse momentum. The angle (ft is the azimuthal angle of k±. The magnitude of 
k± is related to z and y: 



— Pm+l 

= (l-z)y2p a - Ph -k 2 ± (3.7) 
= (1 - z)M 2 e~ t - k\ . 

Thus 

k\ = (1 - z)M 2 e- t . (3.8) 

If we use z and (along with y or t) as our splitting variables, then k 2 ^ is a derived variable. 
Alternatively, we can use k± (and thus k 2 ^ and 0) as splitting variables. Then z is a derived 
variable. 

Lorentz transformation. If y ^ 0, momentum difference p a — p m +i is not exactly equal 
to p a . In order to maintain momentum conservation at each step in the shower, we must, 
therefore, take some momentum from the partons in the final state at the time of the 
splitting. Each parton, with momentum pj, j € {1,2,..., m}, then gets a new momentum 
Pj after the splitting. This includes the momentum pz = pi of the Z-boson. Following the 
shower algorithm of ref. [12] , the momenta pj are determined by a Lorentz transformation, 
pj = Apj with the property 

Pa, + Pb ~ Pm+l = A(p a + p h ) . (3.9) 

However, we use a different Lorentz transformation from that chosen in ref. [12]. 7 If 

p = ap a + (3p h +p± , (3.10) 

then p = Ap is given by 

P = (! + y)«Pa 

Pb (3-11) 



1 

+ 



1 + y 



2p ± k ± k< 

P 7, 1" Ct 



2p a -p h 2p a -p h 



+ P± - ak± . 
An equivalent form that is useful if a 7^ is 

p = (1 + y)ap a 



+ "TTj — \ — o P h ( 3 - 12 J 

(1 +y)a2p a -p h 

+ P± - akx ■ 



7 The choice of Lorentz transformation in ref. [12] takes the needed total momentum from the final 
state partons, but it does not properly absorb the transverse momentum recoil onto the Z-boson. The 
transformation defined here satisfies eq. (3.9), leaves invariant any vector that is orthogonal to p a , Pb, and 
Pm+i and, in addition, transforms pb into a multiple of itself. 
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The momenta p a and py, are not the Lorentz transformed versions of p a and p^. It is, 
however, of interest to know what the Lorentz transformation does to p a and p^. We have 

1 k 2 

Ap a = (l + y)p & + —— ± p h -k ± , 
l+y2p & -p h 

Apb = t~. — Pb ■ 

i + y 

The important feature of this is that the emitted parton m + 1 has transverse momentum 
k± and the momentum Ap a of parton "a" caries the recoil transverse momentum — k±. 
We can understand what happens to this recoil transverse momentum by thinking of the 
shower as proceeding forward in time (oppositely to the way it is generated). Initial state 
parton "a" with momentum Ap a can emit more daughter partons. But a share of the 
transverse momentum — k± is finally transmitted to the Z-boson. 

Let us now look directly at the transformation of the momentum pz of the Z-boson. 
We start with 



Pz = \ e — Pa,+ \ e — Pb+Pz,A 

V S % V S Vb 

= X&PA + X h p B + PZ,± , 



(3.14) 



with 



x h = \ e 



(3.15) 



Then 



Here 



Pz = \ ^ e Y - p a + \ ^ e~ Y - Ph + p Zt± . (3.16) 

V S % V S 17b 

Pz,± = Pz,± ~—k±. (3.17) 
Va 

Thus the Z-boson gets a share x a /r/ a of the recoil transverse momentum. The new rapidity 
is 

r = r + iog(i + ,)--io g (^|A) . (3.18) 

This transformation law is complicated. However it simplifies in the limit that we need 
for this paper. Let us denote by P the Z-boson momentum without its transverse part, 

P = X a pA + X h p B ■ (3.19) 

We are interested in the development of the Z-boson transverse momentum in the region 
p 2 z ± <C M 2 . Therefore we take 

P 2 = M 2 . (3.20) 
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Furthermore, the development of the Z-boson transverse momentum distribution is con- 
trolled by splittings with !/<l. Therefore, we neglect p 2 z ± /M 2 , p 2 ? ± /M 2 , and y compared 
to 1 in Eq. (3.18), giving 

Y = Y (3.21) 
in each splitting. With these approximations, x a and Xb are fixed: 



M 2 Y 
e 



1 M 2 _ y 
[I — e 



(3.22) 



With the approximations eq. (3.22), we have 

Although x a and x\, are fixed, the momentum fractions r\ & and % can change if there is a 
collinear splitting of an initial state parton. 8 

Although we neglect the transverse momentum of the Z-boson in computing its mass, 
we track changes in the transverse momentum as the Z-boson recoils against emissions 
from the initial state partons. For an emission from initial parton "a", the new Z-boson 
transverse momentum is 

X 

Pz,± = Pz,± k± , (3.24) 

as stated in eq. (3.17). 

4. Analysis framework 

The equations of ref. [12] specify quite precisely the evolution of a certain kind of parton 
shower. In order to analyze what the parton shower thus defined produces for the transverse 
momentum distribution of a Z-boson, we develop in this section some theoretical structures 
beyond those presented in ref. [12]. 

4.1 Measurement operators Q 

We are interested in the differential cross section as obtained in the shower at the shower 
final time if, 

In this paper, we will find it useful to represent the desired measurement with the aid 
of an operator Q on the space of statistical states, 

{p±,Y\p(t)) = (l\Q( P± ,Y)\p(t)) . (4.2) 



8 We will see, however, that the real emissions that we need to analyze to follow the development of the 
Z-boson transverse momentum distribution have (1 — z) <C 1. For these emissions, the change in ?7 a and 77b 
is negligible. 
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Here Q(p±,Y) is denned by 

Q(p± , Y) | {P, f, s',c, s, c} m ) 

= 6(pz,±-P±) sfhogf^^) -y) \{ P ,f,s',c',s,c} m ) . ' 

V 2 \pz-paJ ) 

In the subsequent subsections, we will extend this notation in which an operator Q deter- 
mines a measurement on the statistical state. 

4.2 Fourier transformation 

As was noted by Parisi and Petronzio [20], it is useful to analyze the evolution of the 
Fourier transform of the transverse momentum distribution. Thus we study 

/ e- ib ^{p ±) Y\p(t)) = (l\Q(b,Y)\p(t)) , (4.4) 

where 9 

Q(b,Y)\{p, f,s',c',s,c} m ) 

= (2vr)- 2 e-^ S (±1k(?***)-y) \{p, f, s> , d, s, c} m ) . ^ 
\2 \Pz-PaJ ) 

We will refer to (1 Q(b, Y) \p(t)) as the 6-space hadronic cross section. 

4.3 Tracking the momentum fractions and parton flavors 

For our analysis, we will want to keep track of the parton momentum fractions, r; a and r/b, 
and the flavors, a and b, of the incoming partons. Thus we consider the function 

(l|Q(&,y;Tk,77b,a,&)|p(t)) , 

where 

Q(b, Y; 77 a , 77b, a, b)\{p, f, s' , c', s, c} m ) 

= (2^e-^s(llo g ( Pz ^)-Y) (4.6) 
V2 \Pz-PaJ ) 

x 6 a ,aS bt i$(Vz ~ »7a) <%b - Vb) \{p,f,s',c,s,c} m ) ■ 

4.4 Evolution of the perturbative statistical state 

The statistical state vector p(t)) according to our definition contains a factor 

/a/A(% ) M 2 e-*)/ 6/fl (T ?b ,A^e- t ) 



An c (a)n c (b) 2r] ll r] h pA-pB 



(4.7) 



This factor gives the parton-parton luminosity. Here fa/A^la,, P 2 ) and fb/BiVb,^ 2 ) are 
parton distribution functions and n c (a) and n c {b) are the number of colors carried by 



9 We use the same letter, Q, for three different operators, Q(p±,Y), 2(6, Y) defined here, and, in the 
following subsection, Q(b, Y; r) & , 77b, a, b). It should be clear from the context which operator is intended. 
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partons of flavors a and b, namely 3 for quarks and 8 for gluons. We define an alternative 
state vector |p P ert(£)) in which this non-perturbative factor is removed: 

{{p, f, s',c,s,c} m \p(t)) 



/ aM (r/ a ,M 2 e-*)/ 6/B (%,M 2 e -*) , , (4.8) 

= 4n c (a)n c (6)2 WA .p B (W> ^ , S , C } m |p pert (i)) . 

A convenient notation for this is 

=^(i)|p pert (i)) , (4.9) 
where .F(i) multiplies by the parton distribution factor, 

m\{pj,s',c', s ,c} m ) = ^^^y^^^ } • ( 4 - 10 ) 

The evolution equation for |/0 pe rt(O) can b e determined from the evolution equation 
(2.9) for \p(t)). We have 



dt 



\p P ert(t)) +Ht)j t WeAt)) = [H T (t) - V(t)}F(t)\ Ppelt (t)) , (4.11) 



so 



dt 



|p P ert(t)) = ^(t) _1 [Wj(t) " V(i)]jF(t)|p p e rt (t)) 



J-(t)^ 1 



|Ppert(*)) • 



We can write this as 



dt 



\p pc At)) = [n p I cvt (t)-v^ t (t)]\ PpcTt (t)) 



(4.12) 



(4-13) 



Here the revised real and virtual splitting operators are 

nf r \t) = F^Hiit)^) , 

V pcrt (t) = F{t)~ l V(t)F{t) +F{t)~ l 

d 



dt 



(4-14) 



dt 



In the last line here, we have noted that J-{t) commutes with V(i) since V(t) does not 
change momenta or flavors. 

For the perturbative real splitting operator, we have 

{{P, f, S, C, S, C\ m +1 \Hl(t)\{p, /, S, C, S, C\ m ) 

3 (aH(6) /a/A(^a, M 2 e-')/ S/B (r) b , M 2 e-*) 



(4.15) 



n c (a)n c (6) r) a % /a/xfaa, M 2 e -')/ fe/B (7 ?b , M 2 e"*) 

x ({P; /, c', s, 5}^+! |7^f ert (i) /, s', c', s, c}-^) . 

The factors in the second line appear in the definition of the matrix element of Tii in 
ref. [12]. Thus to obtain the corresponding matrix element of 7i pcrt (t), we simply omit 
these factors. We will present detailed formulas for ?^ ert (t) and V pert (t) at the point that 
we need them. 
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4.5 The function to study 

The physics that we want to study is contained in the function defined in section 4.3, in 
which b and Y are measured and, in addition, we measure ry a , a and b. As noted in 
the previous subsection, this function contains a nonperturbative factor that specifies the 
parton luminosity. We remove this factor and study 

(l\Q{b,Y;i] a ,r] h ,a,b)\pp CIt {t)) . 

We will refer to this function as the 6-space partonic cross section. Once we have the 6-space 
partonic cross section, we can obtain the 6-space hadronic cross section (l 2(6,7) p(t)) 
by convolving it with parton distribution functions according to 

(1 \Q(b, Y) \p(t)) = > / drj a drj h ' ' 

^Jo Jo 4n c (a)n c (6)2?7 a 77bPA-PB (4.16) 

X (l|Q(&,Y;77a,77b,a,&)|/9pert(£)) . 

This relation is obtained using eqs. (2.3), (4.5), (4.6), and (4.8). 

Our aim is to study how the 6-space partonic cross section develops as the shower time 
t increases. At shower time 0, it is determined from the Born cross section, as in eq. (2.8), 

(l| Q(&,Y; 77 a , a, &)|ppert(0)) = 12 a Q a bX & x h <5(r/ a - x a ) 5(r] h - x h ) , (4.17) 

with 

[M* y [M* -y <* x 

x a = \ e r , x h = \ e (4.18) 



s 



and 



Q ab = a = gor& = g, 

_ . [l-4|e a |sin 2 (M] 2 + l , Ah , ( 419 ) 

Qab = S 1 . 2 a ^ g and b ± g . 

16sm z (^w) cos z (^w) 

Note that the partonic cross section at t = vanishes unless a is a quark or antiquark flavor 
and b is the corresponding antiflavor. There is no dependence on b because the correspond- 
ing transverse momentum dependent cross section is proportional to a delta function of the 
transverse momentum. As the shower evolves, we expect (l | Q(b, Y; rj^, r]^, a, b) \p pC vt(t)) to 
develop some dependence on b. 

We begin the study of the evolution of this function in the next subsection by outlining 
some key ideas that will go into the derivation. 

5. Outline of the derivation 

We are now in a position to outline the derivation that follows, leaving out most of the 
details and certain subtle points. One of the subtle points is the running of a s . For the 
purposes of this section, we consider a s to be constant. 



- 14 - 



log(fci/M 2 ) 



k\ ~ 1/6 2 




Figure 2: Integration region for initial state gluon emissions. The horizontal axis is the shower 
time; thus emissions are generated from left to right. The vertical axis is log(fc^/M 2 ). The leading 
region for emissions lies between the lines labelled (1 — z) ~ 1 and (1 — z) ~ e~*. The partonic 
6-space cross section for emissions with k\ ^> 1/b 2 approximately vanishes. Thus real gluon 
emissions, indicated by small circles, occur near or below the horizontal line h\ ~ 1/b 2 . 



5.1 Emission kinematics 

We will find that the transverse momentum of the Z-boson is the primarily the result of 
recoils against emission of soft gluons from an initial state quark or antiquark. As we 
saw in section 3, an initial state splitting can be described by splitting variables {t, z, cj)}. 
The shower time t gives the virtuality of the splitting. The momentum fraction of the 
emitted gluon is (1 — z), so that for soft gluon emission we have (1 — z) <C 1. The 
transverse momentum k± of the emitted gluon has azimuthal angle (j) and magnitude given 
by eq. (3.8), k\ = (1 - z)M 2 e - t . 

It is useful to represent splittings by points in the plane of t and log(fc 2 _/M 2 ), as in 
figure 2. An allowed splitting has (1 — z) < 1, so log(fc 2 _/M 2 ) < —t. That is, allowed 
splittings are represented by points below the line labelled "(1 — z) ~ 1" in figure 2. One 
should view this line as having a thickness of order 1. As we will discuss in some detail, 
for a given t, the splitting probability dP from ref. [12] has a term with 

dP~4C F °^dt dz t , (5.1) 
2ir 1 — z + e 1 

including a factor 2 due to having two incoming partons that can radiate. This means that 
splittings are not simply concentrated along the line (1 — z) ~ 1, but are spread over the 
region below this line. In fact, a splitting probability proportional to dtdz/{l — z) would 
give splittings uniformly distributed in log(fc 2 _/M 2 ) and t. However, there is an effective 
cutoff at the line (1 — z) ~ e _t . Thus the splitting probability is approximately constant 
in the region e~* < (1 — z) < 1 indicated in figure 2. 
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We are interested in the partonic 6-space cross section. For this quantity, a real emission 
produces a factor exp(ifcj_ -6), which simply comes from taking the Fourier transform to get 
to 6-space. For this reason, the line k 2 ^_ ~ 1/b 2 in figure 2 is significant. If we integrate the 
real emission probability over an interval of k± in the region k 2 ^ S> 1/b 2 , the factor exp(ifcj_- 
b) averages to zero. That is to say, maintaining a non-zero (l| Q(b, Y; rj a , t%, a, b) p pe rt(i)) 
requires not emitting gluons with k 2 ^ S> 1/b 2 . In the region k"j_ < 1/b 2 , gluons can be 
freely emitted into the final state. We represent gluons that might be emitted in a typical 
event contributing to (l| Q(b, Y; r/ a , r/b, a, b) \p P cvt(t)) by filled circles in figure 2. 

Of course, if there is an approximately uniform probability of emitting gluons in any 
differential unit of area d\og{k 2 _/ M 2 ) dt, then the probability that no gluons are emitted 
for k 2 ^ 3> 1/b 2 is small when b 2 is large. It is approximately 



e-> ^exp^-4C F ^A) , (5.2) 

where A is the area of the triangle in figure 2. This Sudakov factor gives, approximately, the 
6-dependence of the partonic 6-space cross section that we seek. In the following sections, 
we fill in the details of this argument and make it more precise. We will find that the more 
precise analysis leads to a Sudakov factor similar to that in eq. (5.2), but with running 
coupling effects included and an extra "subleading" term. 

5.2 Strategy 

We study the evolution of the 6-space partonic cross section defined in section 4.5. Using 
eq. (4.13), the 6-space partonic cross section evolves according to 

— (l|Q(b,y;r 7a ,r7 b ,a,6)|p per t(i)) ^ 

= (l|Q(b,y;r ?a , % ,a,6)[W I pcrt W-V pcrt (i)]|Ppcrt(t)) ■ 

Our aim is to use suitable approximations to turn this equation into a differential equation 
of the form 

j t (l\Q(b,Y;i h ,r] h ,a,b)\p pcTt (t)) » -K(t,b)(l\Q(b,Y; Vab ,i] h , a,b)\ Ppcrt (t)) . (5.4) 
This differential equation has the solution 
(l | Q(b, Y; 77 a , r/ b , a, b) |p per t(*)) 



exp (- J dr K{t, b) \ (1 1 Q(b, Y; r? a , r/ b , a, b) |pp Crt (0)) 



(5.5) 



Here the initial 6-space partonic cross section has the simple 6-independent form given 
in eq. (4.17). We will see that the evolution of (l 2(6, Y) p(t)) stops for t > t c = 
log (b 2 M 2 e 27E /4), where 7e is the Euler 7. Here t c is approximately the shower time 
at which the lines (1 — z) ~ 1 and k\ ~ 1/& 2 in figure 2 meet; we will see later the reason 
for the adjustment factor e 27E /4. Then 



rtc 

-S 



e~ b = exp I - J dr K{t, 6) J (5.6) 
is the Sudakov factor, for which eq. (5.2) is a simple approximation. 
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5.3 Approximations 

We will need certain approximations to turn eq. (5.3) into the differential equation (5.4). 
We describe these approximations in general terms here. 

First, we note that the behavior of the Z-boson p± distribution for p\ <C M 2 is 
controlled by the 6-space partonic cross section for large b 2 . Thus we are interested in the 
5-space partonic cross section in the region 1/(6 2 M 2 ) <C 1. Therefore, we simply neglect 
l/(b 2 M 2 ) compared to 1 everywhere. 

Second, we neglect e~ t compared to 1. To justify this, imagine letting the system 
evolve from time to a time t\ and then from t\ to t c . Let t\ be large enough so that 
e - * 1 <C 1, but small enough that we can treat t\ as not being a large logarithm. Then 
evolution from to t\ is an approximate version of perturbation theory and gives order a" 
corrections to the Born cross section with no large logarithms. We ignore these corrections. 
For the evolution from t\ to t c , the approximation e~* ~ is justified. Furthermore, we 
can add back the evolution from to t\ using the approximation e~* ~ 0, adding more 
order a™ corrections with no large logarithms. Then we have evolution from to t c with 
the approximation e~* ~ at the cost of changing the result by a" terms with no large 
logarithms. 

For the same reason, we neglect y, eq. (3.4), compared to 1 and k 2 ^_/M 2 compared to 

1. 

Finally, in section 9.2, we will analyze the structure of the splitting function near the 
line k 2 ^ = 1/b 2 . For this analysis, we will make what might be called a low density ap- 
proximation. For initial state emissions, according to eq. (5.1), the density of emission 
points per unit dt and dlog(fc 2 L /M 2 ) is proportional to a s . We treat a s as small. Con- 
sider, then, two emissions, one with parameters {t±, log(fe 2 -,/M 2 )} and the other with 
parameters {t2,log(k 2 l 2 /M 2 )}. Suppose that ti > t\ and that each of these emissions 
is not far from the line k 2 ^ = 1/b 2 . Since the density of points is small, the distance in 
the {t, log(fc 2 /M 2 )} plane between any two points is typically large. This suggests that 
we can neglect e - ^* 2- * 1 ^ compared to 1. To see whether this is justified, suppose that we 
modify the shower algorithm so that it is not allowed to have two splittings be close to the 
line k 2 ^ = 1/b 2 and close to each other. More precisely, choose a distance parameter do 
and require that no two splittings have | log(fc 2 L & 2 )| < do and simultaneously ti — t\ < do. 
Choose do such that e~ d ° is small enough to be neglected but do is not so large that it 
constitutes a large logarithm. Then it is valid to replace e"^ 2- ' 1 ) by zero for any splittings 
that are not excluded. 

The exclusion prescription can be constructed in a different way. Generate 1 — > 2 
splittings without restriction, but add a new splitting function that describes a 1 — ► 3 
splitting: incoming parton "a" emits partons 1 and 2 in the region |log(fc 2 _6 2 )| < do and 
t2~ti < do. The probability for this new splitting should be negative and just big enough 
to cancel the probability for parton "a" to first emit parton 1 and then emit parton 2 in 
this region. One then needs to add similar 1 — ► n + 1 splittings for n > 2 to cancel the 
probabilities to have more than two splittings that are to close the line k\ = 1/b 2 and 
to one another, but we limit this discussion to the simple 1 — > 3 case. The new 1 — > 3 
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splitting probability is proportional to a 2 and the excluded area, cZ§. With our choice of 
do, this area is not large. 

The new splitting function will modify the Sudakov exponent by adding a term pro- 
portional to a 2 and to the length of the line fc^ = 1/b 2 between the two limits in figure 2, 
namely At « log(6 2 M 2 )/2. 

We conclude that replacing iterated splittings by zero in the region in which the ap- 
proximation e - ^ 2 ^' 1 ) ~ is not good results in modifying the Sudakov exponent by terms 
of order a 2 log(6 2 M 2 ). What we actually do is to replace these splittings by what the 
inaccurate approximation e _( -* 2_<1 ^ ~ gives. As long as this approximation results in a 
fractional change of order 1 in the iterated splitting probability, we also modify the Sudakov 
exponent by terms of order a 2 log(6 2 M 2 ). 

The true QCD Sudakov exponent, as discussed in section 12, has an expansion in 
powers of a s . In the coefficient of a 2 , the term with the most powers of log(b 2 M 2 ) is a 
constant times a 2 log 3 (6 2 Af 2 ). The next-to-leading term is a constant times a 2 log 2 (6 2 M 2 ). 
Thus terms of order a 2 log(fo 2 M 2 ) are third-to-leading. The low density approximation 
discussed here changes these terms in an uncontrolled way. 

We have not analyzed here the effect of 1 — > n + 1 splittings for n > 2 that are induced 
by this approximation. However, it should be clear that these induce a™ log(6 2 M 2 ) changes 
in the Sudakov exponent. 

6. Evolution of the partonic cross section 

We can now begin to simplify eq. (5.3), which gives the evolution of the 6-space partonic 
cross section. 

The operator V pert (i) acting on a partonic basis state \ {p, f, s', c', s, c} m ) does not add 
a new parton or change the parton momenta or flavors. For this reason, Q(b, Y; r] a , rj b , a, b) 
commutes with V pert (t). Thus eq. (5.3) becomes 

j t {l\Q(b,Y; Va , Vb ,a,b)\p pcit (t)) = (l\Q{b, Y; rj a , r/ b , a, 6)Wj Crt (t) 

- VP crt (t)Q(b, Y; 7? a , Vb , a, b)\ Ppcit (t)) . 

The operator 7i^ ei (t) does add a new parton and does change the parton momenta. 
Thus it does not commute with Q. To analyze what happens, we break W pert (£) into three 
parts: 




Here Wfs(^) is the part of 7ij CI (t) that generates the splittings of final state particles. 
We have dropped the "pert" notation here because, in the definition of ref. [12], a final 
state splitting does not involve a factor of ratios of parton distributions, so that the part 
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of Ti.f ert (t) that creates a final state splitting is the same as the part of Tti(t) that creates 
the same final state splitting. In the second term in eq. (6.2), we let Hi er \t) be the 
part of 7if cvt (t) that generates the splittings of the incoming parton from hadron A. We 
have decomposed this operator further as an integral and sum of l~L^ cr (t; z, (ft, f), which 
generates the splittings of the incoming parton from hadron A in which the momentum 
fraction and azimuthal angle of the splitting are z and (ft, respectively and the flavor of the 
emitted parton is /'. (For our analysis, /' = g is the most important choice.) Similarly, 
TC^ CTt (t; z, (ft, f) generates the splittings of the incoming parton from hadron B. In a similar 
way, we can divide the virtual splitting operator into three parts, 

V pcrt (t) = VfsW + fdz f vr'^z^j*) 



i 



V (6.3) 



dzj ^£vr(*;^/0 . 

6.1 Final state splittings 

Consider first the effect of a final state splitting. Using the definitions of ref. [12], we find 
that a final state splitting replaces one final state parton by two daughter partons, but its 
effect on the rapidity and transverse momentum of the Z-boson is negligible. When parton 
I with momentum p\ splits to form daughter partons I and m + 1 with momenta pi and 
Pm+i, an amount of momentum Ap = pi+p m+ x —p\ must be taken from the other final state 
partons. The splitting can be characterized by the virtuality variable y = Prp m +i/prQ, 
where Q = p a + pt>- With the splitting kinematics of ref. [12], the needed momentum is 

Ap = -(l-X)p l + (l-\ + y)^-Q , (6.4) 

where 



\=J(l + y) 2 -y%~l-y(-^-l)+--- . (6.5) 
V pi-Q \pi-Q J 

Note that 1 — A oc y for j/ < 1. The momentum Ap is supplied by applying a Lorentz 
transformation to each final state parton i with % ^ {l,m + 1}: = h^vP\- In particular, 
the momentum pz of the Z-boson is transformed with p^ z = A.^ u p u z . 

As discussed in section 5.3, it suffices to consider only splittings with y 1. For y <C 1, 
Ap is proportional to y and hence A^ — 5„ is also proportional to y. Thus the change in 
the Z-boson momentum is small. In particular, the rapidity of the Z-boson changes very 
little, by an amount proportional to y. We can neglect this small change. 

Evidently, the change in the Z-boson transverse momentum must also be small, but 
this statement is not helpful because we are trying to track small changes in the Z-boson 
transverse momentum. To see what happens, we note that the needed transverse momen- 
tum is 

Apj_ = -(l-A)pj,_L • (6.6) 

A fraction of this transverse momentum will come from the Z-boson. Now, we are studying 
the evolution of the probability that the Z-boson transverse momentum pz ± is small and 
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(6.9) 



remains small. For this to happen, the transverse part of pi must be small, of order pz,± 
or smaller. 10 Thus, recalling that (1 — A) is of order y, we have 

Pz,± ~Pz,± ~ Cyp z ,± , (6.7) 

where C is of order 1 or smaller. That is, the fractional change in the Z-boson transverse 
momentum due to a final state splitting is negligible. 

This discussion can be summarized by saying that, to a sufficient approximation, 

Q{b, Y-^rjk, a, b)H FS {t) ~ Hps^QfrY;^,^ a, b) . (6.8) 

This gives 

— (1 1 Q(b, Y;r] a , r] h , a, b)\p pcrt (t)) 

~ (l|[^Fs(0 _ VFs(t)]Q(b,Y;r]a,rib,a, b)\ 

Ppert 

(*)) 

+ fdz [ n (i|Q(&,F;»7a,%,a,6)wr rt (i;^^/0 

- Vr\t, z, </>, f)Q(b, Y; 7? a , 77b, a, b)\p pelt (t)) 

Jo J -it Z7T 

- V b pert (t; z, 0, f')Q(b, Y; 77a, 77b, a, b)\p pcrt (t)) . 

This is useful because the definition of Vfs(^) is based on the requirement, designed to 
insure that the shower conserves probabilities, that 

(l|[Hps(t)-VPs(t)]=0 • (6.10) 

Thus the first term in eq. (6.9) vanishes. We must analyze initial state splittings, but we 
can ignore final state splittings entirely. 

6.2 Initial state splittings 

We now turn to initial state splittings. We relate Q applied after the splitting to Q applied 
before the splitting. The relation is 

Q(b, Y; 77a, 7? b , a, b)Hi ert (t; z, <j>, f) « Hi eTt (t; z, 0, /') exp(i ^ b kj_j 

x zQ(b,Y;zf) a ,rib,a- f',b) . 

To understand this, first look at the b dependence, using the definition (4.6) of Q. When 
we apply the operator Q after the splitting, it produces a factor exp(— ibpz j_) where 
pz j_ is the Z-boson transverse momentum after the splitting, which is related to the Z- 
boson transverse momentum after the splitting, pz,±, and the transverse momentum in 
the splitting, k±, by eq. (3.24), 

Pz,± = Pz,± ~—k±. (6.12) 

7?a 



We discuss this in section 9.1. 



- 20 - 



Recall that k± is specified by z and <fi: it has azimuthal angle 4> and square k"j_ = (1 — 
z)M 2 e~ t . Thus 

exp(— ib-pz±) = exp(— ibpz ±) exp ( i — b k± ) . (6.13) 

V 'b: J 

The factor exp(— ib-pz,±) is generated by the Q operator before the splitting but the 
second factor in eq. (6.13) must be supplied. The dependence on the momentum fractions 
is simple. According to eq. (3.6), the momentum fraction r?b is unchanged by the splitting, 
while he momentum fraction ?) a after the splitting is related to the momentum fraction ry a 
before the splitting by rj a ~ zr/ a . Thus the r/ a argument of Q before the splitting is zf] & and 
there is a jacobian factor z because Q is defined with a delta function. Finally, the flavor 
a after the splitting is related to the flavor a before the splitting by a = a — /', where we 
use the notation u — g = u, g — u = u, etc. 

With these observations, our equation for the variation 6-space partonic cross section 
with shower time is 

^ (l| 2(6, Yifjz, fjt,, a, b) |p pert (t)) 
/* j'^ d(j) 

~ J o dz J 2tt^2 ( 1 |^a(i;^,0,/';b,^;?7a,'%,S,fe)|Ppert(t)) ^ ^ 

+ dz —^(l\lC h (t]z,(j)j']b,Y]fj a ,fj h ,a,b)\p pcrt (t)) , 

JO J-TT Z7T ^7 

where /C a describes a splitting of an initial state parton "a" from hadron A and is given by 
(l|/C a (t; z, 4>, /';&, I"; r/ a , r? b , a, b)\p pcrt (t)) 

= {l\Hr\t; z, <P, /') expfi ^ 6-fc ± ) Z Q(6, Y; zfj a , a - /', 6) (6.15) 

- Vr\t; z, <f>, f')Q(b, Y; 7? a , 5, 6)|p P ert(i)) • 

The operator /Cb f° r a splitting of the initial state parton from hadron B is the same with 
the roles of "a" and "b" are interchanged. 



7. The real splitting function 

At this point, we need to know some details about the operation of "H a crt (£; z, <f>) on a 
general partonic state Up, f, s', c', s, c} m ) ■ Fortunately, we need only the inclusive splitting 
probability 

(l\Hr\t-zAJ')\{pJ,s\^s,c} m ) . (7.1) 

Formulas from refs. [12, 13, 14] for this quantity are reviewed in appendix A. The most 
important case to consider is that of a q — > q + g splitting orag-> q + g splitting. However, 
we include all flavor choices. We treat separately two kinematic regimes: (1 - 2) < 1 and 
(1 — z) ~ 1 since the results in these two regimes have rather different structures. Our 
subsequent analysis of evolution will make use of the results in these two regions. 
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The operator 7i\ ert contains a factor a s . In ref. [12], the argument of a s was denoted 
by fi R and left unspecified. In general, n R can be a function fjt R (z, t) of the kinematic 
variables that describe the splitting. Our default choice in this paper is 

^ 2 R = X R (l-z + y)M 2 e- t , (7.2) 

where 

/ Ca[67-3tt 2 ] - lOraA 

>[- [ 3(33-2„ t ) J ' < 7 ' 3 ' 

Except when (1 — z) is very small, this is approximately the constant \ R times k 2 ^ = 
(1 — z)M 2 e~ t . In section 12, we will see why choosing [i R approximately proportional 
to k 2 ^ is useful and we will see why the choice given in eq. (7.3) for the proportionality 
constant is also useful. 

Leaving these points for later, we can immediately understand why a factor (1 — z + y) 
in eq. (7.2) is preferable to the simpler choice (1 — z). As we will see, having a running 
scale fi R as in eq. (7.2) with either a factor (1 — z + y) or a factor a factor (1 — z) affects 
the Sudakov exponent that we obtain in section 12. When a s ((x R ) is expanded in powers 
of a s (M 2 ) one obtains terms proportional to logarithms of (1 — z + y) or (1 — z) times 
extra powers of a s . These terms improve the matching between the Sudakov exponent 
obtained with the shower and the true QCD Sudakov exponent. We will find that, with 
the accuracy of matching that we can obtain, logarithms of (1 — z + y) or of (1 — z) are 
equivalent. However, we can still ask which is more desirable in general. The running 
a s (n R ) incorporates some features of the singularity structure of higher order graphs into 
the leading order splitting. We note that the leading order splitting kernel from ref. [12] 
has a singularity 1/(1 — z + y). This suggests that the higher order contributions might 
naturally contain a logarithm of this same variable, {1 — z + y). In contrast, the leading 
order splitting kernel does not have a singularity when (1 — z) — > at fixed y, so it would 
not be natural to introduce a logarithm of (1 — z) into the expansion of a s (fi R ) . Indeed, soft 
gluon emissions correspond to y — > and (1 — z) — > together, while (1 — z) — > at fixed 
y corresponds to an anticollinear emission in which incoming parton "a" emits a gluon in 
the direction of incoming parton "b". There are such singularities, but they are associated 
with emissions from incoming parton "b" rather than from incoming parton "a" . For this 
reason, our default choice (7.2) for [i R has a factor (1 — z + y) rather than (1 — z). 

7.1 (1-z) < 1 

We use the splitting operator for an initial state splitting from ref. [12], using the splitting 
variables t, z, (ft defined in section 3 of this paper. When the emitted parton is a gluon, the 
splitting probability has a "soft gluon emission" singularity that corresponds to a factor 
1/(1 — z) when y <C (1 — z) <C 1. In this section, we extract the terms that have this soft 
gluon factor; other terms will be included in the following subsection, where we study the 
regime (1 — z) ~ 1. We note, in particular, that contributions in which the emitted parton 
is not a gluon do not give a 1/(1 — z) contribution. 
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We use the results of ref. [12], particularly eqs. (12.20), (12.21), and (12.22), as reviewed 
in appendix A. 11 Part of eq. (12.21) is a function An-. There is some freedom in choosing 
this function. We use the definition in eqs. (7.2) and (7.12) of ref. [14], which are equivalent 
to eqs. (A. 6) and (13.4); we discuss the reason for this choice in section 13.2. We neglect 
y compared to 1 and (1 — z) compared to 1. However, we do not neglect y compared to 
(1 — z). In the shower formalism of ref. [12], when parton "a" splits by emitting a gluon, 
we include interference graphs. The soft gluon is emitted from parton "a" in the amplitude 
and by some other parton k in the complex conjugate amplitude, or by parton "a" in the 
complex conjugate amplitude and by parton k in the amplitude. Here parton k could be 
parton "b" or could be a final state parton. For this reason, there is a sum over the indices 
k of helper partons, with k ^ a. 

After some calculation, we find in appendix A that 



{l\zUr\t-zA- f'))\{ P ,f,s',c',s,c} m ) 

~ S f',S J2(i S '}rn\{s}m)(-l)({c'} m \T k -n\{c} r 



(7.4) 



2ir 1 — z + y 



f{z,y,<t>;r k ) ■ 



Here is the rapidity of the helper parton k relative to the rest frame of p a + and 
is its azimuthal angle; then / is 



\ ~ 1 r k 2 \/( 1 - Z )V (i , n . 2r k 2 V 

f{z,y,<pr k )~ : cos(# - fa) + e k - — 

1 — z + y 1 - z + y 



(7.5) 



Because we are calculating an inclusive splitting probability, indicated by the measurement 
function (l|, there is a quantum inner product ({s'} m |{s} m ) between the spin states in 
the amplitude and the conjugate amplitude. 12 There is also a color inner product, which 
is non-trivial because it contains a color matrix T£ for emitting a gluon from line "a" and 
a color matrix T£ for emitting the gluon from line k, summed over the eight colors c of 
the emitted gluon. This same factor appears in the Catani-Seymour dipole subtraction 
formalism for next-to-leading order perturbative calculations [21]. 

This splitting function is complicated because of the function /. To understand /, 
denote the rapidity of the emitted gluon relative to the rest frame of p a + by r. For a 
soft gluon, (1 — z) -C 1 and y <C 1, we have 

(1 - z) « e 2r y . (7.6) 

For a given splitting time t, y is fixed and we integrate over z. There is a near singularity 
in this integral for (1 — z) — > 0, but this near singularity is cut off when (1 — z) becomes 
comparable to y. That is, r > 1 is favored in the integration. Writing e 2r y for (1 — z) in 

n We note that the first factor on the right hand side of eq. (12.22) of ref. [12] needs to be complex 
conjugated; the * is missing. 

12 In a spin averaged shower, {s'} m = {s} m at every step and this factor is 1. 
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f(z, y, (f>\ r k ) gives 



f{z,y,<f>;r k ) 



1 + e zr 1 + e zr 



(7.7) 



We see that when r ~ r&, all three terms in f(z, y, <p; r k ) are comparable, so that we have 
quite a complicated function. However, 

f(z, V, & rk) ~ 1 when r»r fc . (7.8) 

We thus see that when the rapidity of the emitted gluon is large compared to the 
rapidities of previously emitted gluons, the splitting function simplifies to 

[l\zHr t {t-,z,cl>\f)\{pJ,s',c!,s,c} m ) 

« <V,g E <{^ m |M4(-l)<{cV|T fc -T a |{c} m > . (7 - 9) 

fc^a 

Now the only dependence on the helper parton index k is through the color factor. This 
enables us to perform the color sum as described in ref. [12], 

-E({ c 'M Tfe ' Ta K c W = <{c'} m |T a -T a |{c} m > = C a ({c'} m \{c} m ) , , T1Q , 



where 



\C F a^g 

Ca={ ^ ■ (7.11) 

\C A a = g 



Thus 

(l\zHl eit (t;z,<f>; f')\{ P J, s',c',s,c} m ) 



&f* (W}m\{s} m )({c'} m \{c} m ) C * 



That is, using eq. (2.5), 



2-7T 1 — z + y 



«b(a4) 2 , lr , , \ (7-13) 

^f'^^TT^Ca- —(l\{ P ,f,s',c',S,c} m ) . 



;i|^f rt (t;z,0;/')|{p,/, S , ,c / , S ,c} m ) 

2tt 1 — z + y 



7.2 (1 - z) ~ 1 

Now consider the collinear limit (1 — z) ~ 1. We use the general result from ref. [12], 
reviewed in appendix A. We neglect e - ' compared to 1 as discussed in section 5.3. As 
described in appendix A, the result is a simple color structure that multiplies the standard 
(unregulated) DGLAP splitting kernels P a ,a'{z) and the ratio of the number of colors for 
parton flavor a' = a + / to the number of colors for parton flavor a, 

{l\zHl e *{t-,zAJ')\{p,f,s'^,s,c} m ) 



n, 



:(« + /') «sK) 1 p , W1|/ . , , , v ( 7 - 14 ) 
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8. The virtual splitting function 



We now look at the virtual splitting function. We need 

{l\Vl^(t-z,cj>j')\{p,f,s',d,s,c} m ) . 

Eq. (4.14) gives 



J.l 



(l|V pert (t) = (l| 



v{t) + Htr 1 Jt Ht) 



1.2) 



Following ref. [12], we define the virtual splitting operator V(i) from the requirement 

(l|V(t) = , (8.3) 

which guarantees that shower evolution does not change the cross section when we integrate 
over all final states that start from the hard scattering. We break H\ into pieces as in 
eq. (6.2), 

Hi(t) = H FS (t) + jf dz f [Wa(t; z, 4>, /') + H h (t; z, <t>, /')] • (8.4) 

This differs from eq. (6.2) in only one respect: the parton distribution factors that were 
omitted in the pieces of Wj >ert (t) are included in the pieces oiHi(t), as indicated in eq. (4.15); 
for a splitting of parton "a", using fj a « r] a /z and a = a + /', the relation is 



(l\Ha,(t; Z, <j), f')\{p, f, s', c, s, c} rn ) 

_ f(g+f')/A{Va./z,M 

n c (a + f')f a/A ( Va ,Mh 



{l\zHl ert (t;z,<t>J')\{p,f, s',c',s,c} m ) 



(8.5) 



We can similarly write V(i) as a sum and integral in the form 

V(t) = Vps(t) + J\z f i^Y, z, 0, /') + Vb(t; z, 0, /' 

With this notation, the definition in ref. [12] is 

(l|V a (*; = (l|Wa(<; • 
The needed matrix element involving the derivative of J-(t) is 



(8.6) 



(8.7) 



(ll^t)- 1 -^)!^, /, s',c',s,c} t 



ifa/A^AFe-') |/6/ B (r7a,M 2 e-*) 



+ 



/ aM (77a,M3 e -*) f b/B ( Va , M^) 

x (l|{p,/,s',c' ,s,c} m ) . 



- 25 - 



The parton distribution functions obey the DGLAP equation 



L3 



dt 



fa/A{Va,M Z e 



2-t\ 



Pa,a+f'{z) f {a+ /')/a(W Z, M 



2„-t\ 



5f 



f',s 



2C a 
1 - z 



la 



f a/A {r ] ^M*e- t )\+0{a 



(8.9) 



Here P a ,a{ z ) are the standard (unregularized) DGLAP kernels, C a is either Cf or Ca as in 
eq. (7.11), and 



la 



1 Cf a ^ g 



^[llC A -2n f ] a = g 
Following our default choice (7.2) for the argument of a s , we have used 



.10) 



=a s (A R (l-z + y)M 2 e"*) 



in the evolution equation for the parton distribution functions. This is not the standard 
choice, but it can be accommodated without changing the parton distribution functions 
by modifying the evolution kernel at next-to-leading order and beyond. That is, the terms 
indicated by 0(a 2 s ) are modified from what they would have been had we used a s (M 2 e _i ). 
These terms do not affect our analysis. 

With these results, we can wr ite the complete V pert (t) as a sum and integral in the 
form used in eq. (6.3), 

v^\t) = v|f (t) + fdz [ v a pcrt (t; z, 0, /') + v b pcrt (t; *, 0, /')] . (s.ii) 



For Va° r (t; z, 0, /')) we have 

{l\Vr\t;z,<P,f')\{p,f,s',c',s,c} m ) 



nc(a) f( a +f')/A{vJz,M 
n c (a + Df a/A (v a ,M^ 



^ (1 \znr\t; z, cj>, f')\{p, f, s', c', s, c} m ) 



x (l|{p,/,s / ,c' > s,c} m ) . 
This simplifies in two limits, which we now discuss. 



2Cg 

1 - z 



la 



.12) 



13 If the parton distributions used in shower generation obey the second or higher order DGLAP equation, 
then there are more terms with extra powers of a s . We do not display possible higher order terms in order 
to keep the notation from becoming complicated, but in the end we will find that their inclusion would 
change the Z-boson transverse momentum distribution at a level that is beyond the accuracy that we aim 
for. 
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8.1 (1 -z) <C 1 

Consider the limit (1 — z) <C 1. Then the leading terms in Va° rt contain a factor that equals 
1/(1 — z) as long as 1 — z y; contributions without this factor can be neglected. There 
is no such factor unless /' = g. Thus only /' = g is important. In the factors with ratios 
of parton distributions, a + f = a and ry a / 'z r/ a ; thus these factors are well approximated 
by 1. This gives 

(i|vr t (t;^0,/o|{p./y.^*»c} ro ) 

« 5 fjS (1 1 ^P ert (t; z, 0, /') I {p, /, s', c', s, c} m ) 



2tt 



j-Pa,a(2) - Y~, +7a|(l|{p,/,S / ,C / ,S,c} m ) . 



U3) 



14) 



In the second line on the right hand side of this equation, P a ,a{z) contains a term pro- 
portional to 1/(1 — z). Other terms, which do not contain this factor, can be ignored. 
Similarly, we can ignore the term — j a since it has no 1/(1 — z). In fact, the term in P a ,a( z ) 
proportional to 1/(1 — z) is 2C a /(l — z), which cancels the term — 2C a /(l — z). This leaves 
the very simple result, 

{l\vr n {t V zAJ')\{pJ,s'^,s,c} m ) 

«*/'.g(l| ^r rt (t; z, <j>, /') I {p, f, s',c', s, c} m ) 
8.2 (1 - z) ~ 1 

Now consider the collinear limit (1 — z) ~ 1. For the matrix element of 7i\ ert in this limit, 
we use eq. (7.14). This gives 

(l\Vi eit (t;z^j')\{p,f,s',c',s,c} m ) 
_ f(a+f)/A(ri a /z,M 2 e- t ) « s (^r) 



4/A(la,M 2 r' 

«s(/4) 



2vr 



Pa,a+J>{z) 



f(a+f)/A{r] a /z,M 2 e 



2Cg 

1 - z 



la 



.15) 



x (l|{p,/, s',c',s,c} m ) . 
The two terms involving ratios of parton distributions cancel, leaving the very simple result, 
(l\vr\t-z^J>)\{pJ,s\d,s,c} m ) 



8ft 



a 



Ml 

2tt 



2C a 
1 - z 



la 



{l\{p,f, s',c',s,c} m ) 



.16) 



9. Evolution for t < t c 



The lines k\b 2 ~ 1 and (1 — z) ~ 1 in figure 2 cross at a shower evolution time given 
approximately by t = log(b 2 M 2 ). At this time, the nature of the evolution of the 6-space 
partonic cross section changes. We define a critical shower evolution time t c by 



t c = log Qb 2 M 2 e 27E ^ 



(9.1) 
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Here 7e is the Euler constant. We will see later the reason for the factor exp(27E)/4. In 
this section, we make use of the results from the previous section to analyze the evolution 
before t c , e~* 3> e~ tc . 

Recall from eq. (6.14) that that the evolution is specified by operators that we called 
/C a and /Cb, 

~ (l| Q(b, F;?7a, 77b, a, b) \p pei t(t)) 



dt 



I dz [ W ^V(l|/C a (t;z>,/ / ;fo,y;7 ?a ,r ?b ,a,6)| Ppcrt (t)) (9.2) 

J J -IT /7T fl 



f 

+ (b splitting term) . 

The operator /C a is defined in eq. (6.15). We will analyze /C a ; the analysis of /Cb is the 
same. 

In the end, /C a is very simple. However, we will have to analyze it in stages, pruning 
away complications at each stage. 

9.1 Distribution of final state partons 

We are concerned with the evolution of the 6-space partonic cross section, 

(l\Q(b,Y;r] a ,rj h ,a,b)\p pcrt (t)) . 

This is an inclusive quantity, but we first note something about the structure of the partonic 
final states that contribute to it. Recall from its definition (4.6) that Q contains a factor 

e ~ 11 e ' (9.3) 

i>i 

where Pj ; ± is the transverse momentum of the jth final state parton (other than the Z- 
boson). When we form the fr-space partonic cross section, we integrate over all of these 
transverse momenta. The integration region in which one or more partons have p 2 ± 3> 1 jb 2 
gives a negligible contribution because the exponential factor exp(ipj j_ • b) averages to zero. 
It is as if Q contained a factor Y\j0(p 2 ± < C/b 2 ), where C is a constant that is large 
compared to 1 but with logC not large. That is, Q effectively projects onto partonic final 
states in which no partons have been emitted in the region above the line k 2 ^ ~ 1/b 2 in 
figure 2. 

One simple consequence of this concerns any real parton splitting that has the potential 
to change the flavor a or momentum fraction r) & of parton "a" . Such a splitting must be 
close to the line (1 — z) ~ 1 in figure 2. To see this, note first that in order to significantly 
change r/ a it is necessary that (1 — z) not be tiny compared to 1. Second, in order to change 
the flavor a, the parton emitted into the final state must not be a gluon, but the splitting 
functions for these splittings do not have (1 — z) — ► singularities and hence have negligible 
probabilities of occurring with (1 - z) < 1. 

We now note from eqs. (3.8) and (9.1) that 

k\ = 4e 2 ^ l -^e t <- t . (9.4) 
b z 
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Thus any emission with (1 — z) ~ 1 and e~ t S> e~' c has k 2 ^ ^> 1/b 2 . We conclude that all 
emissions that contribute to the 6-space partonic cross section for t < t c leave the flavor a 
of parton "a" unchanged and leave its momentum fraction rj a approximately unchanged. 
Thus the initial conditions 

??a = x & , 77b = x h (9.5) 

remain (approximately) true and the flavors a and b of the incoming partons do not change. 
With eq. (9.5), eq. (3.23) for the splitting variable y reads 

y = e- t . (9.6) 

In the following section, we note another consequence of the fact that Q effectively 
projects onto partonic final states in which no partons have been emitted with k 2 ^ S> 1/b 2 . 

9.2 Angular ordering 

In the range e~ l S> e~* c and k 2 ^ ~ 1/b 2 , both real and virtual emissions contribute to the 
evolution. For e~ t 3> e~ tc and k 2 ^ ^> 1/b 2 , only virtual emissions contribute. In either 
case, it appears that the analysis of the evolution will be very complicated. The virtual 
emission operator Vi ev (t; z, 4>, g) is given in eq. (8.12) in terms of the real emission operator 
H.& CTt (t; z, cp,g). The real emission operator is given in eqs. (7.4), (7.13), and (7.14). Here 
is where the complications are. There is a certain algebraic complexity and there is a 
non-trivial color structure. Worse, there is a sum over helper partons with labels k, with a 
different splitting function for each k. Thus, in order follow the evolution of the inclusive 
quantity (b, Y\pP eit (t)), it appears that we need to track the structure of the complete 
partonic state. 1 

Fortunately, there is a simplification available. We consider the real or virtual emission 
of a gluon with k\_ > 1/b 2 at shower time t with e - * S> e~ tc . The rapidity r in the Z-boson 
rest frame of the gluon is approximately given by 

*+^og(SVlog(|r) (9-7) 



2 to V M2 . 

as long as e~ l -C 1. Using eq. (9.5), this becomes 



Curves of constant r are shown in Figure 3. 

The rapidity of a previously emitted gluon (that is, a real gluon in the final state) is 

1 (k 2 \ 

r k = t k + 2 \og ^ . (9.9) 



14 This would not be quite so bad in the leading color approximation generally used for parton showers. 
Then, we would need only the momentum of the parton that is color connected to the incoming quark or 
antiquark, which would generally be the gluon previously emitted from that incoming parton line. 
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Figure 3: Integration region for initial state gluon emissions, as in Figure 2. The slightly curved 
lines indicate lines of constant rapidity for the emitted gluon, with rapidity near at the left and 
large rapidity on the right. 

Here we have omitted the log (1/zk) because any previously emitted gluon would have had 
(1 — Zk) -C 1, so log(l/zfc) ~ 0. Since the gluon was previously emitted, we have tk < t. 
Since it was a real emission, the transverse momentum satisfied k 2 ^ k ~ 1/b 2 . On the other 
hand, k 2 ^ > 1/b 2 . Additionally, log(l/z) > 0. From these inequalities, we conclude that 

r k < r . (9.10) 

If either k\ » b 2 or k\ k <C b 2 , we have 

e 2r >e 2r '= . (9.11) 

In the event that k\ and k 2 ^ , are both of order 1/b 2 , this strong inequality may not hold. 
However, the emission probability per unit t and per unit log(k 2 _/M 2 ) is small, of order 
a s . Thus we can apply the low density approximation developed in section 5.3 to conclude 
that in this event we can approximate 

e 2< > e 2tk (9.12) 

at the cost of affecting terms in the Sudakov exponent that we will obtain at the level of 
third-to-leading terms and beyond. Eq. (9.12) then implies eq. (9.11). 

For this reason, we can make the approximation eq. (9.11) in the expression (7.4) for 
Hr\t;z,<P,g). This gives the much simpler formula eq. (7.13), for which the details of 
the state of the shower at time t are not needed. 

We are now prepared to write an evolution equation for (l | Q(b, Y; r] a , a, b) \p peT t(t)) 
for t < log(5 2 M 2 ). We consider first z <C 1, then z ~ 1. Then we combine these two cases. 
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9.3 Contribution from (1 - z) < 1 

In this subsection, we consider fC a in the region 



l-z«l , e -t > e" tc . (9.13) 

That is, we consider splittings that are represented by points well below the line (1 — z) ~ 1 
in figure 2 and well to the left of t = t c that passes through the lower right vertex of the 
triangle. In this region, gluon emission from line "a" is important because the splitting 
function has a term with a 1 /(l — z) factor. However, splittings with f'^g lack this factor 
and may be omitted. With this approximation, eq. (6.15) becomes 

(l|/C a (i;z, (f), /'; b, Y; rj & , rj h , a, b)\p pcrt (t)) 



3 /',g 



1 \Ur\t; z, 0, g) exp ( i ^ bk±) zQ(b, Y; z Va , 7? b , a, b) (9.14) 

V zr ]a, J 

- Vr\t- z, 0, g)Q(b, Y; Va , rtb, a, b)\p pcrt (t)) . 



Since 1 — z <S 1, we can also replace z by 1 in the exponential factor and in the argument 
of Q(b, Y; zr\ a , 7/t,, a, b). We can also use eq. (9.5) to replace x a /rj a in the exponent by 1. 
Finally, we use eq. (8.14) for Vr rt - With these chan ges, we have 

(l\JC a (t; z, 4>, f; b, Y; 7/ a , rj h , o, b) |p pe rt(*)) ab ^„ 

■ ^ g (9.15) 
« 8 r>g (1 \zHr\t; z, </>, g)Q(b, Y; r/ a , Vh , a, b)\p pcrt (t)) [exp(i bk ± ) - 1] . 

We can make one more simplification. The right hand side of eq. (9.15) vanishes when 
k\ -C 1/fa 2 . In the alternative case that k\ > 1/b 2 , we can apply the lesson of section 9.2 
so as to use the simple form of (ll^a^ given in eq. (7.13), 

(i|^^(t ; ^;g)|{p,/yy, Sl c} m ) 

~^4^C a - 1 —- i (l\{p,f,^,s,c} m ) . 

2tt 1 — z + e _t 

Here we have replaced y by e~ t according to eq. (9.6). Thus 

(l |/C a (i; z, <j), /'; b, Y;r] a , i] h , a, b)\p peit (t)) 

wtf 1 [exp(i6.fe ± )-l] (9.17) 

2ir 1 — z + e 1 



x 



(l | Q(b,Y; Va,Vb, a, b)\p per t(t)) . 



Of course, the expression that we have used for (l|zWa Crt is not accurate in the region 
k^_ -C 1/b 2 ; nevertheless, both the exact form and the new form vanish in this region, so 
this inaccuracy is not a problem. We can use eq. (9.17) for both k 2 ^ <C 1/b 2 and k 2 ^ > 1/b 2 
so long as (1 — z) -C 1 and e~* ^> e~ tc . 

Notice the key feature of eq. (9.17) that the object for which we seek an evolution 
equation, 

(l\Q(b,Y;r] a ,rj h ,a,b)\p peit (t)) , 



- 31 - 



appears on the right hand side of this result. 

We can make one more simplification in this. Recall that we need the average over the 
emission angle <fi of (l j/C a |/o pe rt(*)) • The only (f) dependence is in the factor exp(i bk±). 
When we take the average over (f>, we get a Bessel function: 



2vr 



l|£ a (i; z, 0, /'; b, Y;r/ a , rj h , a, b)\p peit (t)) 



2ir 1 — z + e 1 

(l|Q(&, Y;r] a ,r] h ,a,b)\p peTt (t)) . 



We will want to integrate this over z. Consider the integral between limits z\ and 
Z2, with e~ t <C (1 — z\) < (1 — z-x) <C 1. The denominator (1 — z + e - *) can then be 
approximated by simply (1 — z). Similarly, the argument of a s is, using eq. (9.6) and then 
e-* < (1 - z), 

fjfc = A R (1 -z + e~ t )M 2 e- t fa A R (1 - z)M 2 e~ t = \ R k 2 ± . (9.19) 

Using = (1 — z)M 2 e~ t , we can change the integration variable to k\. This gives 
f z ' 2 f n d<p , 

dZ 2TT ( 1 F a ( t ' 2 '^' /''' & >^>' r ?a,?7b,a,&)|Ppert(*)) 

-wj'f ^wi^ii-D-!] <»»> 



X 



'l| Q(&, F; 7? a , 77b, a, b)\p pcrt (t)) 



Thus we encounter an integral of the form 

rQl fju2 

F(QlQi)= -r^f(k 2 ± )[J Q (\k ± \\b\)-l] , (9.21) 

Jq'{ k ± 

where the function f{k\) allows for the dependence of a s on k\. As we have previously 
argued, if Q 2 and Q 2 are both much smaller than 1/& 2 , then we can replace [Jo(|fcj_| |6|) — 1] 
by zero since Jo(|fcj_||6|) ~ 1 for |fcj_||6| — ► 0. Similarly, we have argued that if Q 2 and 
Q\ are both much larger than 1/& 2 , then we can replace [Jo(|fcj_||6|) — 1] by —1 because 
•A)(|fc±||fr|) is a rapidly oscillating function in the integration range. But what happens if 
Q 2 <C 1/b 2 while Q 2 S> 1/b 2 ? There is a simple approximation that is accurate provided 
that f{k 2 A ) is a slowly varying function. We can replace [Jo(|fcj_||6|) — 1] by — > 
2e _7B ), where is the Euler constant. Thus the integral becomes 

rQl Hk 2 

F(Ql Q 2 ) appr0 x = - / %t /( fc D • ( 9 - 22 ) 

The integrals in eq. (9.21) and (9.22) differ by terms proportional to a power of l/(Q 2 b 2 ) 
or a power of Q 2 b 2 and by a term proportional to the second derivative of f(kjj with 
respect to log(fc 2 L ) in the region near k"j_ ~ 1/b 2 . The exact statement can be found in 
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appendix B. Note that the second derivative of a s (AR,fcj_) with respect to log(fc^) is of 
order o^Ar^) 3 , so that this is a good approximation for our purposes. 
With this approximation, 

(l|/C a (t; z, 4>, /'; b, Y; r/ a , rfc, a, b)\p pcr t{t)) 

~ ~ Sf* ^ C a r -^ Zi 6{kW > 4e" 2 -) (9-23) 

X (l|S(6,y;7/ a ,7/b,0,6)|/>pert(*)) , 

where the ~ in this case indicates that this approximation works inside the integration 
over 4> and z. 

Specifically, the derivation so far covers a range of z between limits e - * <C (1 — z\) < 
(1 — Z2) <C 1. We have in mind that (I — Z2) is chosen small enough that the approximation 
(1— z) <C 1 applies within the integration range. The region (1 — z) ~ 1 is treated separately 
in the following section. 

Similarly, we choose (1 — zi) = Cie - *, where C\ 3> 1 but log C\ is not large. Then the 
integration range < (1 — z) < C\e~ l needs a separate treatment, which we now provide. 
For that range, the denominator 1 — z + e - ' in eq. (9.18) provides a lower cutoff on (1 — z) 
at around (1 — z) ~ er l '. Thus the integration range that needs a separate treatment is 
really e _i < (1 — z) < C\e~ l . What happens here depends on the factor [Jo(|fej_||6|) — 1] in 
eq. (9.18). Since k\_ = (1 — z)M 2 e~', define the value (1 — Zb) of (1 — z) that corresponds 
to = 1 by 

(1 - z fe )M 2 6 2 e~* = 1 . (9.24) 

If (1 — Z{,) e~ l , then |fcj_||b| 3> 1 in the integration range e~ l < (1 — z) < C\e~ l and 
we can approximate [Jo(|fej_||&|) — 1] — > —1. If (1 — 0&) S> Cie~ t , then |fcj_||6| <C 1 in the 
integration range e~* < (1 — z) < C\e~ l and we can approximate [Jo(|fcj_||6|) — 1] — > 0. 
Both cases are covered by the approximation 

[^o(|fc±||b|)-l]^^(fcib 2 >4e- 27B ) . (9.25) 

In the remaining case, e~* < (1 — zj,) ^ Cie - *, we will use this same replacement. This 
is not an accurate approximation. However, this inaccuracy occurs only near the point at 
the intersection of the line (1 — z) = e~* and the line k\ = 1/6 2 in figure 2. Following 
an argument like that in section 5.3, we recognize that this inaccuracy does not matter 
because it does not lead to contributions with a large logarithm in the Sudakov exponent. 
Thus eq. (9.23) is also a sufficient approximation when applied inside an integration over 
z that includes (1 — z) — > 0. 

9.4 Contribution from (1 — z) ~ 1 

The approximations in the previous section cover the integration region (1 — z) -C 1 for 
e^ t S> e~* c . For the integration region (1 — z) ~ 1 with eT l S> e~* c , only virtual splittings 
contribute: 

(l|X: a (t;2:,^,/ / ;6,F;?7 a ,r/b,a,6)|pp er t(t)) 

« -(l|Vr t (*;^^/')Q(6,l r ^a,»7b,a,6)|ft» rt(t)) • 
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For the matrix element of V a er , we can use eq. (8.16). This gives 



(1 |£ a (i;z, 4>, /'; b, Y; r/ a , r} h , a, b)\p pcrt (t)) 



5f 



Ml 

2vr 



2C a 
1 - z 



la 



[l\Q(b, Y; rj a , rj h ,a, b)\p pcrt {t)) 



).27) 



As in the previous section, we notice that the object for which we seek an evolution equa- 
tion, (l| Q(b, Y; ?7 a , r/b, a, b)\p pcr t(t)) , appears on the right hand side of this result. 



9.5 Evolution for all z 

Compare eq. (9.23), which applies for (1 — z) <C 1, with eq. (9.27), which applies for 
(1 — z) ~ 1. Notice that the theta function 9(k 2 _b 2 > Ae~ 2lE ) is always 1 when (1 — z) ~ 1 
and e _t 3> e~' c , that (1 — z + e _t ) is the same as (1 — z) when (1 — z) ~ 1 and e _i 1, 
and that the term — 7 a with no 1/(1 — z) factor is negligible (compared to 1/(1 — z)) for 
(1 - z) < 1. Thus these two forms match and we can combine them in the form 



'l|/C a (*; z, 4>, f; b, Y; rj a , r] h , a, b)\p pcvt (t)) 



5 f/ 



/',g 



2tt 



2C a 



l-z + e 



7a 



9{k\b 2 >4e~ 27B ) 



.28) 



^l|Q(b,y;r? a ,?7b,a,6)|ppert(i)) • 



This gives the evolution equation, for e '>e ic and assuming that parton "b" caries the 
opposite flavor from parton "a" , b = a, 



d 
dt 



1 1 Q(b, Y; 77 a , ?7 b , a, 6) |p pcr t(*)) 



2vr 



2C a 



1 - z + e" 



7a 



0(fci& 2 > 4e- 27E ) 



.29) 



^1 1 Q(b,F; 77a, r?b, a, 6) |/9 P ort(i)) • 



There is a factor 2 here because emissions from both initial state lines contribute equally 
when b = a. 

This analysis could apply when a = b = g, which would be relevant for computing 
the transverse momentum distribution in Higgs boson production. However, the case of 
interest in this paper is that a is a quark flavor and b is the corresponding antiquark flavor, 
or vice versa. In this case, C a = Cf and "y a = (3/2) Cf. Then 



d_ 



(l | Q(b, Y; 77 a , 77b, a, b) |p per t(*)) a 



dz 



(4) 



2tt 



2C F 



2 

1 - z + e~* 



).30) 



[l\Q(b, Y; 7/a, 77b, a, b)\p pa . t (t)) 
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10. Evolution for t > t c 



After shower time t c , the evolution changes character. To see what happens, it is simplest 
to consider (l | Q(b, Y) \p(t)) instead of (l| Q(b, Y; r/ a , 77b, a, b) Lo per t(i))- We have 

|(i|e(b,y)|p(t)) = (i|Q(6,y)[wi(t)-v(t)]|p(t)) (101) 
« (i|Q(6,y)Hi(t)-v(t)Q(b,y)| P (t)) . 

Consider the relation of Q(6, Y")Wi(t) with 7ii(t)Q(b, Y). As we have noted in our previous 
analysis, each emission from an initial state line with transverse momentum k± produces 
a phase factor exp(i6-fcj_) in the result of applying Q(b,Y) after the splitting; this phase 
factor is not present in the result of applying Q(b, Y) before the splitting. However, all 
emissions for e~* <C e~ tc have k 2 ^ <C 1/b 2 . Thus this phase factor is simply 1. Thus 

Q(b,Y)Hi(t) Kt Hi(t)Q(b,Y) . (10.2) 

This result gives 

j t (l\Q(b,Y)\p(t)) « (l|[Wi(t)-V(t)]Q(6,y)|p(t)) 
= . 

Thus (l Q(6, y)Lo(t)) stops evolving at t = t c . The b-space partonic cross section, 
(l \Q(b, Y; 77 a , r/b, a, 6) |/0 per t(*)) j does evolve for t > t c because initial state emissions can be 
collinear and thus change the momentum fractions and flavors ?7 a ,r7b,a, 6. What happens 
is that the partonic cross section evolves and the parton distribution functions evolve with 
opposite evolution kernels, so that the net change of (l\Q(b,Y)\p(t)) with t vanishes. 

11. Solution at t — t c 

We can solve eq. (9.30) for the 6-space partonic cross section with initial condition (4.17), 
evolving from shower time t = to shower time t = t c . This gives 



(l|Q(6,y;r? a ,77b,a,6)|ppert(^c)) = 12 a Q ab x a x b <5(r/ a - x a ) S(rj b - x h ) 

xexp(-5 (M 2 ,b 2 )) , 

where the Sudakov exponent So is 

2 3 



(11.1) 



S (M 2 ,b 2 ) = [ tC dt Cdz ^^-2C F 
Jo Jo 27T 



l-z + e-t 2 

It is useful to change variables from t to k\ = (1 - z)M 2 e~\ giving 



6(k 2 ± b 2 > 4e- 2 " lE ) . (11.2) 



M 2 ju2 



dk l [\„ a ^i) 



S (M 2 ,b 2 )= / ^ / dz -^ 2 c F 



2(1 -z) 



;i - z ) 2 + k 2 jM 2 2 



(11.3) 



We can now approximate So(M 2 ,b 2 ) in a fashion that makes it simpler. We are 
interested in the behavior of Sq(M 2 , b 2 ) for M 2 b 2 S> 1. In particular, we are interested in 
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terms with logarithms log(M 2 6 2 ). The region of the k\ integration that dominates for large 
b 2 is k 2 ^ <C M 2 . In contrast, the region k 2 ^ ~ M 2 cannot contribute logarithms log(M 2 fo 2 ) 
and is thus not of interest for us. For this reason, we can make the approximation k 2 ^ -C M 2 
inside So(M 2 ,b 2 ). Additionally, it is useful to expand ct^/x^) in powers of a s (fc 2 L ). Then, 
after performing the z-integral, we obtain contributions to the integrand of the k\ integral 
of the form 

a s (klr[log(M 2 /kl)] m . 

The leading terms in this counting are those with m = n. We need to keep track of those, 
along with the next-to-leading terms with m = n — 1. However, we will not be interested 
in the coefficients of terms with m < n — 2. For this reason, we drop all such terms in the 
approximations below. 

To proceed with this program, we note that the renormalization scale in a s is, using 
the definition (7.2) and then eq. (9.6), 

k 2 



Mr 



Ar(1 - z + e- t )M l e~ t = \^k\ \ 1 + 



;i - z) 2 m 2 



(11.4) 



Expanding a s (/i 2 l ) in powers of a s (k^J and displaying only the order a s and a 2 terms, we 



have 



Os(A*r) _ a s {k 2 ± ) 



2tt 



2tt 

-2ft log! 14 



2/3 1 log(A R ) 
k 2 



a s (k 2 A 
2tt 



1 



) 2 M 2 



a s (k 2 
2vr 



+ 0{al 



Here 



33 - 2 m 
12 

The z-integral that multiplies the order a s term is 



dz 



2(1 



(1 - z) 2 + k 2 /M 2 2 



log 



M 2 
k 2 



- +log 1 + 



M 2 



(11.5) 



(11.6) 



(11.7) 



The most important term here is the logarithm log (M 2 /fc 2 L ). The term —3/2 is next-to- 
leading, so we keep it. The third term in eq. (11.7) contributes to the k 2 ^ integral only for 
k 2 ^ ~ M 2 but is power suppressed in the dominant region k 2 ^ -C M 2 . Consequently, we 
neglect this term. This same integral multiplies the order a 2 term proportional to log(Aft). 
Here, we keep the large logarithm log (M 2 /kY) but neglect the term —3/2 in which the 
extra power of a s is not multiplied by a large logarithm. For the remaining a 2 term, we 
need the integral 



dz 



2(1 -z) 



(1- z) 2 + k 2 /M 2 



log 1 + 



(1 



) 2 M 2 



TT' 



+ o 



/fci 

\M 2 



Here there is no logarithm, log (M 2 /k 2 ), multiplying the extra power of a B , so we neglect 
this contribution entirely. An analogous analysis shows that we can neglect entirely the 
order and higher order terms in the expansion eq. (11.5). 
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With these approximations, So ~ S where 



S{M 2 ,b 2 ^ 



dk 2 

4e- 2 ~<E /b 2 



k 2 



2tt 



2C F 



log 




4/3iC F log(A R ) log, 2 



M 2 
1 



(11.9) 



12. Result 

We have seen that the 6-space hadronic cross section [l Q(6,K) Uf)) evolves from t = 
to t = t c = log (b 2 M 2 e 27E /4) and then stops evolving. Therefore the Fourier trans- 
form of the physical pj_-space cross section is given by (l \Q(b, Y) \p(t c )^j. Furthermore, 
(l Q(b,y) jj(t c )) is given as a convolution of parton distributions and the 6-space partonic 
cross section (l| Q(b, Y; r? a , r] h , a, b) |p pe rt(*c)) by eq. (4.16), 



{l\Q(b,Y)\p(t c )) = [ dri, fdm, 



L/A^Ae-^/b^h/d^Ae-^/b 
An c {a)n c {b) 2rj ai i] h p A -p B 



;i2.i) 



x (l|Q(6, Y;r7 a ,T7b,a,&)|ppert(£c)) • 
We have found that the 6-space partonic cross section at t = t c is given by 

(l | Q(b, Y; a, b)\p pcrt (t c )) = !2aQ ab 5(l -x a /r] a )5(l - x h /i] h ) 

x exp(-S(M 2 ,b 2 )) , 



(12.2) 



where the Sudakov exponent S given in eq. (11.9). This has been obtained with the 
approximations discussed in section 5.3 and in section 11. In particular, terms in S of the 
form a™ log n_1 (M 2 6 2 ) are affected by the approximations. 

We can compare this to the result of ref. [17] for QCD. In ref. [17], there are two 
arbitrary parameters, C\ and C2, that do not affect the result summed to all orders of 
perturbation theory, but do affect the perturbative expansion. Here, we take the simplest 
choice, C\ = 2e _7E and C% = \. Additionally, the result is given for an arbitrary choice of 
the scale p 2 = M 2 e~ t that appears in the parton distribution functions. Here t should be 
close to t c . We choose t = t c . With these choices, the result of ref. [17] is that the hadronic 
5-space cross section is given by eq. (12.1) with 



[l\Q{b, Y;ri a ,r] h ,a,b)\p{t c )) 



QCD 



12a Y" ]Q a >b> C a i a [ — ,a s 




4 e -27E 



C&'i 



4e- 2 TE/6 2 k"_ 



A(a s (k 2 ± ))log 



M 2 
k 2 



(12.3) 



+ B(a s (k 2 ± )) 



-37- 



Here A, B, and C have perturbative expansions in powers of ot s of the indicated arguments. 
The first terms in these expansions are [17, 22, 23, 24, 25] 



A(a s ) = 2C F ^+2C F K(^) 2 + 
Ztt \Ztt/ 



4°^ + 
2vr 



B(a s ) -- 
C a > a (z,a s ) = 8 a > a5(l - z) 



197 34n f 20vr 2 8n f vr 2 8C(3) 



+ 



+ 



27 + 3 



(12.4) 



+ 



2^ 



<*a'a ^ 3 (1 - Z) + 5(1 - Z) 



3 2vr 2 23 
2 + ~3~ " ~~ 



6 



+ 5 ag z(l - z) 
Here we follow the notation of [7] in defining 

K = C A 



67 
18 



21 



5 rif 
~9~ 



(12.5) 



We begin the comparison of the shower result (12.2) with the QCD result (12.3) by 
noting that they agree at the Born level, order a®. This agreement was built in when we 
chose the Born cross section as the starting point for shower evolution. 

We next note that the QCD cross section exponentiates in 6-space. This is a statement 
about (l\Q(b,Y; r] a , rjb, a,b)\p(t c )) q CT) . We consider the logarithm 15 of this function and 
expand it in powers of a s (M 2 ) and log(M 2 6 2 ). To obtain this expansion, we need to 
expand a s (/c 2 L ) in powers of a s (M 2 ) and log(M 2 /fc 2 L ) using the evolution equation for 
a s , then perform the integration over in the Sudakov exponent S. Similarly, in the 
coefficients C, we need to expand a s (4e~ 27B /b 2 ) in powers of a s (M 2 ) and log(M 2 /6 2 ). 
The general form of the perturbative expansion would be 



log [(l| Q(b, Y; 77a, r/ b ,a, &)|p(tc)) QGD ] 

a s (M 2 ^ 



£ £ D nm (Y) 



n=l m=0 



2tt 



)5 



(12.6) 



We know that the maximum power m of log(M 2 b 2 ) must be no larger than 2n because 
each loop can give at most two logarithms, one from a collinear singularity and one from 
a soft singularity. If we look at the perturbative expansion of the partonic fo-space cross 
section instead of its logarithm, then terms with m = 2n actually occur. By the statement 
that the shower cross section exponentiates in 6-space we mean that, for all n, 



D nm = for m > n + 1 . (12.7) 

15 In taking the logarithm, we consider the functions C to be operators on the space of parton distribution 
functions. Thus the first term, S a ' a S(l — z), represents the unit operator. The coefficients D nm (Y) in 
eq. (12.6) are then, in general, operators on the space of parton distribution functions. One could diagonalize 
these operators by taking moments, but there is little reason to do so. 
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This property of the QCD result (12.3) is shared by the shower result (12.2). For the 
shower, it is a simple consequence of having a differential equation in which the derivative 
with respect to shower time of the partonic 6-space cross section is a kernel times this same 
partonic 6-space cross section, where the kernel has one a s and one logarithm. 

One sometimes finds discussions of summations of large perturbative logarithms ex- 
pressed in terms of the directly observed cross section, which can be written as a convolution 
of parton distribution functions and a partonic pj_-space cross section defined analogously 
to the partonic 6-space cross section. Thus one may write 



[l| Q(P±, Y \ »7a, ^a, &)|/3pert(*c)) 



n=l m=l 

One can then discuss the coefficients E nm . However, this is not very useful. For instance, 
no condition on the E nm for m = 2n and m = 2n — 1 will imply eq. (12.7). It is for this 
reason that we compare the shower result and the QCD result with respect to the logarithm 
of the partonic 6-space cross section. 

Let us now compare the n = 1 terms in eq. (12.6) between the QCD result (12.3) and 
the shower result (12.2). We first note that the terms D\q do not match between the two 
results: the order a s contribution to C a / a (z, a s ) in the QCD result is lacking in the shower 
result. These terms arise in part from the one loop virtual graphs, which are not included 
in the shower splitting functions. For that reason, the shower does not get these terms 
right. In our analysis we have have, accordingly, completely neglected contributions to the 
evolution kernel that have a factor a s with no logarithms log(M 2 fo 2 ). 

We now compare the al terms with one or two powers of log(M 2 fo 2 ). For the shower 
result (12.2), we have 

4 

D12 = ~ , 

(12.9) 

D u = - — [log2-7E] -4 . 

This is the same as D\2 and Du in the QCD result. It is not surprising that the leading 
coefficient, D\2, matches. The contribution D\2 log 2 (M 2 6 2 ) is 4Cf times the area of the 
triangle in figure 2. The contribution Du log(M 2 b 2 ) is more subtle. This contribution is 
associated with the edges of the triangle. The bottom edge of the triangle is associated 
with the lower limit on the k 2 ^ integral, which is 4e _27£; /b 2 instead of simply 1/6 2 . In 
the shower evolution, this value arises as a property in integrals involving Bessel functions 
Jo(|fc_!_||f>|). In the analysis of the contribution to evolution associated with fc^ near this 
limit, it was important that the rapidity of a potential gluon emission was large compared 
to the rapidities of previous gluon emissions. In the analysis of the contribution associated 
with the left hand edge of the triangle, it was important that the denominator of the 
splitting function is 1 — z + e _t and not, say, 1 — z + 2e - '. This denominator is related 
to the treatment of the soft gluon interference diagrams. Finally, the right hand edge of 
the triangle is associated with the term —3/2 in eq. (11.9) for the Sudakov exponent. This 
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term comes from the — 7 a in the virtual part of the DGLAP evolution equation for the 
parton distribution functions, eq. (8.9). 

Let us now look beyond the order a s contribution to the Sudakov exponent. 

We note that the QCD result and the shower result share the contribution inside the 
k 2 ^ integral, 

In the shower result, this form reflects our choice to use defined in eq. (7.2) as the 
argument of a s ; this is a choice that is not required by the logic of the shower, in which 
we factor softer interactions from harder interactions. With this choice, we match the 
contribution in eq. (12.10) between the shower result and the QCD result. Expanding in 
powers of a s (M 2 ), we see that 

D%T = D mn D n>l, m = n + l . (12.11) 

We note also that the QCD result and the shower result share the contribution inside 
the kj_ integral, 

V ±J 2C F — . (12.12) 



2vr \ 2, 

This term contributes to the coefficients D nm for n = m. However, the QCD result has a 
contribution 

(a s (k 2 ± )\\ (M'< 



2C F jnog -2- . ( 12 - 13 ) 



which is to be compared to 




2tt / °U 



M 

4/3 1 C F log(A R ) log( p-) . (12.14) 



It is not a surprise that these do not match for Ar = 1, since we have used a leading 
order shower and this is a second order contribution. However, this term in the Sudakov 
exponent contributes to the coefficients D nm for n = m and n > 2. If we would like to 
match these terms, we can, by means of a trick introduced in ref. [7]. We choose Ar such 
that 

-4/3iC F log(A R ) = 2C F K . (12.15) 

This is what we have done in our default choice given in eq. (7.3). For rif = 5, this is 
Ar ~ 0.41. If we take this default choice, then 

<n Wer = ^n D n>l, m = n . (12.16) 

This scale choice applied to q — > q + g splittings produces the second order term in 
A(a s ) for Z-boson production. The same formalism applies to the transverse momenta 
of Higgs bosons produced in hadron collisions. In that case, the incoming partons are 
two gluons instead of a quark and an antiquark. The Sudakov exponent has the same 
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form, but with different functions A, B, and C. For A, the first two terms are now 
2CA(a s /(2vr)) + 2C J 4i^(a s /(27r)) 2 , but the ratio of the first and second coefficients in A(a s ) 
is unchanged [4] . Thus the same factor Ar applied to g — > g + g splittings produces the 
second order term in A(a s ) for the Higgs boson transverse momentum distribution. Of 
course, this is a trick that applies for the particular purpose at hand. It is not equivalent 
to using a second order splitting kernel. 

For the Higgs boson p± distribution, we also replace B\ = — 2^ q = —4 by B\ = — 27 g , 
where -y g is given in eq. (8.10) [4]. 

13. Other choices 

In this section, we briefly investigate the Z-boson transverse momentum distribution that 
would result if we were to make other choices for the construction of the shower. 

13.1 Spin-averaged, leading-color shower 

We have organized this paper as an analysis of the shower evolution equation in ref. [12], 
with some minor modifications. This evolution equation contains the effects of spin corre- 
lations and applies to arbitrary color states. It contains the physics that we believe should 
be approximately represented in a parton shower event generator. However, the nature of 
the evolution is such that finding a way to represent the equations in a practical computer 
program is challenging. To start with a base approximation, one can average over the spins 
and take the standard leading color approximation. We analyzed this case in ref. [13] . With 
a spin average and a leading color approximation for each splitting, the evolution equa- 
tions have the right form to be implemented as a Markov process, as is commonly used 
for parton showers. 16 We are thus led to ask whether the results of this paper still hold if 
one averages over spins and takes the leading color approximation at each parton splitting 
step. The answer to this question is that the results do still hold. 

Spin does not matter because the observable of interest, the Z-boson transverse mo- 
mentum, is independent of parton spins. The spin dependence is represented in the spin 
factor in eq. (7.4), ({s'} m |{s} m ): the parton spins {s} m in the quantum amplitude must 
equal those in the quantum conjugate amplitude, {s'} m , but the splitting probability is 
independent of what the spin values are. If we average over spins, this factor goes away. 

Color is of some importance. However if we use the exact color dependence then 
eq. (7.10) allowed us to reduce the color dependence to a simple factor of Cf or Ca- If 
we use the leading color approximation, then each gluon is treated as carrying color 3 3 
instead of color 8. The allowed parton pairs forming a dipole are then the color connected 
partners. Then eq. (7.10) still holds with Cf replaced by Ca/2. One can get back to 
Cf by simply adjusting the quark-quark-gluon couplings in the splitting functions. With 
this understanding, the results of the previous section calculated with full color remain 
unchanged in the leading-color approximation. 

16 We are, in fact, currently engaged in writing computer code for this purpose. 
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13.2 Alternative choice of dipole partitioning function 

Within the framework of a virtuality ordered shower as presented in ref. [12], there is a 
choice to make. Consider the emission of a soft gluon from a color dipole composed of 
the initial state parton "a" and one other parton, k. The total emission probability is 
proportional to 

H ak ({p} m +l) + HU{P} m+ l) = - • (13.1) 

V" Pm+l 'Pa Pm+1 'Pk 

Here we use the eikonal approximation, which is valid for very small p m +i- The shower 
algorithm divides this into pieces, H &k and H ka , defined by 

Hak{{p}m+i) = — A' ak ({p} m+1 ) - , 

vr Pm+rPa, Pm+VPk ^ ^ 

H k a{{P}m+l) = — A' ksL ({p} m+1 ) — . 

7T p m +l -Pa, Pm+l 'Pk 

Here A' &k and A' ka are both positive and obey 

Ai k + A' ka = l . (13.3) 

The contribution proportional to A' k is considered to be an emission from parton a and 
is treated using the momentum mapping associated with emission from parton a. The 
contribution proportional to A' k is considered to be an emission from parton k and is 
treated using the momentum mapping associated with emission from parton k. 17 Thus the 
function A' ak tells how to partition the soft emission from a dipole into two parton splitting 
contributions. 

In this paper, we use the A' &k defined in eq. (7.12) of ref. [14], 

A l = Pm+l-Pk Pa-jPa+Ph) ^ ^ 

ak Pm+l-Pk Pa' (Pa+Pb) + Pm+l 'Pa Pk '(Pa +Pb) 

with A' k& given by the same expression with a <-> k. In the p a + pb rest frame, this is 

A > = l-cosfl m+ i, fc 

afc (l-cos0 m+life ) + (l-cos0 m+ i, a ) 

We have seen that with this choice the shower will produce the proper structure for the 
Drell-Yan transverse momentum distribution. 

One can consider using instead the A' &k defined in eq. (7.13) of ref. [14], 

A ' ak = -„ p ? +1 '? k . (13.6) 

Pm+l-Pk + Pm+l' Pa. 

This is the form used by Catani and Seymour and by others for the purpose of partitioning 
the dipole splitting probability in defining perturbative dipole subtractions. The most 



17 In eq. (12.21) of ref. [12], the splitting probabilities are expressed in terms of functions A^k and Aka 
which are related to A'^ k and A' ka by eq. (7.2) of ref. [14]. 
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important difference between the two choices is that A' ak as given by eq. (13.4) is invariant 
under the scaling pk — > Xpk- Thus this A' ak depends on the angle between the emitted 
gluon and the momentum of parton k, but not on the energy of parton k. 

We have investigated the possibility of using A' k defined by eq. (13.6) to generate 
the parton shower. We find that if we do so, the argument given in sections 7.1 and 9.2 
fails. Now the probability of emitting a gluon m + 1 from the initial state parton "a" with 
virtuality smaller than that of a previously emitted gluon k can depend on the momentum 
Pk- As a result, we no longer obtain a differential equation for the 6-space partonic cross 
section that has only a kernel times the 6-space partonic cross section on the right hand 
side. The Z-boson will get less transverse momentum recoil from emitted gluons, but we 
cannot say how much less. 

We can draw a lesson from this: a parton shower built on splitting functions that reflect 
the collinear and soft singularities of QCD does not automatically sum large logarithms 
that appear in a physical cross section of interest. Seemingly minor details like the choice 
of A' ak matter. 

13.3 Catani- Seymour dipole shower 

The shower discussed in this paper is the virtuality ordered shower of ref. [12] with an 
improved momentum mapping for initial state radiation and with the renormalization 
scale for a s given in eq. (7.2). A close relative is the Catani-Seymour (CS) dipole shower, 
proposed in ref. [26] and implemented in refs. [27, 28]. Here one takes the splitting functions 
and momentum mappings of the Catani-Seymour algorithm [21] for defining subtractions 
in NLO perturbative calculations and uses them instead to define a shower algorithm. We 
examine this choice here, with the evolution variable taken to be 



as in refs. [26, 27, 28]. 

The most important difference between the shower of this paper and the CS dipole 
shower lies in its treatment of momentum conservation. Let us consider the (backward) 
evolution of an initial state parton by emission of a gluon into the final state. The initial 
state parton is on shell with zero transverse momentum. It is replaced by a prior on- 
shell parton with zero transverse momentum and the final state gluon that has non-zero 
transverse momentum. Thus there is a transverse momentum imbalance. In the shower 
of this paper, the recoil from the final state gluon is taken by all of the other final state 
partons collectively, by means of a Lorentz transformation k% = Aki. With this choice, we 
have seen that the recoil transverse momentum is predominately taken by the Z-boson. 
In the CS shower, for each splitting there is a partner parton, the other parton in the CS 
dipole. Let us suppose that we use the leading color approximation, which is commonly 
applied in parton showers. Then the partner of a radiating initial state quark is the parton 
that is color connected to the quark. For the first gluon emission, the partner of incoming 
parton "a" is the other incoming parton "b". In this case, the CS scheme assigns the 
recoil transverse momentum to the Z-boson. However, once one gluon has been emitted, 




(13.7) 
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the partner of "a" is the previously emitted gluon. Then the CS scheme assigns the recoil 
transverse momentum to the partner parton. That is, the previously emitted gluon recoils 
against the newly emitted gluon and the Z-boson does not take any recoil. 

This leads to a curious situation. With a transverse momentum ordered shower, the 
shower first evolves by not emitting any gluons. The probability not to emit any gluons 
up to a certain shower time t± is a Sudakov form factor. At some point, one real gluon is 
emitted. Then evolution of the Z-boson transverse momentum distribution stops because 
the recoil from further gluons emitted is taken by one of the gluons previously emitted. 
This problem has been analyzed and addressed in ref. [29]. 

We can formulate the issue in the style of our analysis of the previous sections if we 
supply one additional variable in the measurement operator Q, a variable recoil that takes 
values T and F. Including this variable, our operator expressing possible measurements is 
Q(b, Y; rja,, 7/b, a, b; recoil). If recoil = T, then when a gluon is emitted from either of 
the initial state partons, the Z-boson will share in the recoil against the gluon transverse 
momentum. If recoil = F, then, when a gluon is emitted from an initial state parton, the 
Z-boson does not share in the recoil. 

The evolution equation of the recoil = T cross section is 

-^(l\Q(b,Y; Va ,rj h ,a,b;T)\ Ppcrt (t ± )) 

di±_ (13.8) 
= -(l|VP firt (<)Q(b ) y;» 7a ,i7 b> a J 6;T)| i0 pert(t±)) • 

This does not involve the recoil = F cross section. Using the approximations developed 
in the preceding sections, we can solve this equation. We find 



[l|e(6,y;^,o,6; T)\ Ppcrt {t ± )) 

= 12 a Q ab x a x b <5(?7a - x a ) 5(r] h - x h ) 



2 dp 2 a s (X R p 2 ) 
x exp - / — ^- - 2C F 



p 2 2vr 



log 



~W)~2 



(13.9) 



Note that for very large t±, the recoil = T cross section tends to zero. 

The evolution equation of the recoil = F cross section has only the recoil = T cross 
section on the right-hand side: 



dt 



'l|Q(6, Y;va., Vb, «> b; F)\p peTt (t±)) 



[ iz : £i?""( I ;*-* ± ) 



/' 

x {l\zHr\t- z, <P, /') Q(b, Y; z^, Vh , a - /', 6; T) |p pert (t ± )) 
+ (b splitting term) . 



[13.10) 



The total partonic 6-space cross section is 

Q(b, Y; r/ a , r? b , a, b) = Q(b, Y; r/ a , Vh , a, b; T) + Q{b, Y; Va , r, h , a, b; F) . (13.11) 



- 44 - 



For very large t±, the recoil = T contribution vanishes and we are left with the recoil = 
F contribution, which we can obtain by solving eq. (13.10) with the use of eq. (13.9). The 
result for large t± is 

(l | Q(b, Y;r] a , rj h , a, b; F)\p pcrt (oo)) 



Jo rl 



x 12 a Q a -f>,b Xa%b S(zrj a - x a ) 5(r? b - x b ) 



x exp 



— 2 2G F 



fe 2 /r 2vr 
n c {o) a s (/4) 1 



log 



7? 



(13.12) 



n c(a - f) 2tt z 
+ (b splitting term) . 



P ( y\ 2zC F fy, g 2zC F fy, g 

^a-f',a{Z) ; h 



1 



1-z + e- 



We see that the partonic cross section does not exponentiate in 6-space. In fact, it is sim- 
plest in pj_-space instead: if we take a Fourier transform back to p± , the factor Jo(|fcj_||&|/;z) 
provides a delta function 2(27r) 2 <5(p 2 L — k^/z 2 ). This leaves a relatively simple expression 
consisting of a Sudakov factor and a splitting function for the single allowed splitting. 

Again, we see that a parton shower built on splitting functions that reflect the collinear 
and soft singularities of QCD does not automatically sum large logarithms that appear in 
a physical cross section of interest. Seemingly minor details like the momentum mapping 
matter. 



13.4 Transverse momentum ordered shower 

The shower examined in this paper is based on virtuality ordering. What would happen 
if we kept everything else the same but used transverse momentum ordering? That is, 
what if we replace the evolution time t of eq. (3.3) by t± defined in eq. (13.7)? We keep 
the momentum mapping the same as in the virtuality ordered shower, so that the transfer 
of recoil transverse momentum to the Z-boson is the same as in the main body of this 
paper. Transverse momentum ordering is used in Pythia [30, 31], but the algorithm that 
we consider in this section differs in some respects from that of Pythia. 

With k± ordering, the first stage of shower evolution is simple. For t± <C log(M 2 6 2 ), 
only virtual splittings contribute to the evolution of the 6-space partonic cross section; 
the rapidly oscillating Jo(|fc±||b|) factor that multiplies TCj Crt eliminates the real emission 
contribution. 

The later stages of shower evolution are also simple. The argument of section 10 implies 
that the rate of change of the 6-space hadronic cross section vanishes for t± ^> log(M 2 b 2 ). 

We are left with the region t± ~ log(Af 2 b 2 ). This is the bottom of the triangle in 
figure 2. With a virtuality ordered shower, this region is sampled from left to right, which 
is in order of increasing rapidity of emitted gluons - that is, decreasing emission angles. 
Then the splitting function in eq. (7.4) took the simple form in which the complicated 
function / was replaced by 1 because the rapidity of the real or virtual gluon was much 
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larger than the rapidity of previously emitted gluons. With t± as the evolution variable, 
there is no guarantee that a previously emitted gluon had smaller rapidity than a new 
real or virtual daughter gluon. Thus the form of the evolution equation changes and the 
evolution of the partonic 6-space cross section depends on the previous emission history. 

Because the evolution is more complicated than we have dealt with, we do not know 
what the result is. We do note that the difficulty arises only for a limited region, t± ~ 
log(M 2 6 2 ). 

13.5 Angle ordered shower 

What would happen if we used angular ordering instead of virtuality ordering? That is, 
what if we replace the evolution time t of eq. (3.3) by tz = — log(tan(#/2)) where 6 is the 
splitting angle. Then ordering in tz is ordering in splitting angle, starting from large angles 
and progressing to smaller angles as tz increases. At the same time, we would use a revised 
version of the dipole partitioning function A' lk , as specified below. This gives us two of the 
main features of Herwig [32, 33], although if we start with the shower described in this 
paper and change only these features, we certainly do not have exactly Herwig. 

To be more precise, let P = x a p a + x^pb be the Z-boson momentum approximated as 
having no transverse part but as obeying P 2 = M§. Consider the splitting of a parton 
with label I and momentum pi and let n\ be the lightlike vector 

Ml 

m=Pi-^pP- (13-13) 
Let the daughter partons have momenta pi and p m +i- Then define 

"-Nf ^p'" ) ■ < 1314 » 

2 \p m+ i-pi P-niJ 

For an emission from initial state line "a" , tz is the rapidity of the splitting that we used 
in section 9.2. From eq. (9.8) we have, for y <C 1 and r/ a ~ x a , 

^ =t+ ^ iog (S) +i ° g G) • (i3 - i5) 

Curves of constant tz are shown in figure 3. Over most of the region of the figure, (1— z) <C 1 
so log(l/z) ~ 0. This gives straight lines in the t-log(fc 2 _/M 2 ) plane. As these lines 
approach the line (1 — z) ~ 1 in the figure, the log(l/z) term becomes important and the 
tz curves bend. 

The curves of constant tz extend down to arbitrarily small virtuality, M 2 e~ t . In 
a virtuality ordered shower, one would simply stop the evolution when M 2 e~ t reaches 
some infrared cutoff at which application of a hadronization model is more appropriate 
than continued use of perturbative showering. For the angle ordered shower, we should 
impose this same cutoff as a upper limit on integrations over t at fixed tz- Equivalently, 
if we parameterize the curves of fixed tz using k"j_, then we impose a lower limit on the 
integrations over k 2 ,. 
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One should do more than simply switching the evolution variable. One can also modify 
the dipole partitioning function A' , as used in eq. (13.2). For an angle ordered shower, it 
is standard [6] to begin with 



A 



(0)/ 



1 



afc 2 

A (0)> = I 
fca 2 



1 COS fl m+ l, a ~ COS 9 m +l,k 
1 - COS a ,fc 

_ cos0 m+ i )a - cos6* m+ i ifc 



(13.16) 



1 - COS a ,fc 

where the angles are defined, for instance, in the Z-boson rest frame. These functions add 
to one, but are not everywhere positive. One can then replace the functions -ff a ,fc and -fffc.a 
defined in eq. (13.2) by their averages over the azimuthal angles defined by rotating p m +i 
about, respectively, the directions of p a and As shown in ref. [6], this gives a simple 
result, 

v \ ° s /i(o)Vr-i \ 1 l-cos0 a ,fc 

7T dK E m+1 I 1 ~ COS 0m+l,a) (1 ~ COS 9 m +l,k) 

a s 1 0(0 a ,m+l < ^a,fe) 
^ VT l-COS0 m+ l, a ' 

tr ^ " s /i(°)Vr»-i \ 1 l-cos0 a ,fc 

7T Ka E^ +l (1 - COS m+ l, a ) (1 - COS m+ l,fc) 

«s 1 0{0k,m+l < 0a,fc) 

7T £^ +1 1 - cos m+ i,fc 
The procedure is equivalent to replacing A^? 1 and A^' by 



^afe — #(#a,m+l < ^a.fc, 
A ka, = 0{Qk,m+l < Oa,k. 



1 — COS 0- 



1 - COS a ,fc 
1 - COS m+ i, a 
1 - COS a ,fe 



(13.18) 



These functions are positive. They do not sum to 1, but, as just noted, one gets the right 
result after suitable averaging over azimuthal angles. 

It is now easy to see how the partonic 6-space cross section Q(b, Y; rjg,, 77b, a, b) evolves 
as tz increases. At each step dtz, there are contributions from all values of fc^ down to the 
infrared cutoff. We divide the integration over into three distinct regions: k 2 ^ 3> 1/b 2 , 
k 2 ^ ~ 1/fo 2 , and k 2 ^ <C 1/b 2 . In the region k\ -C 1/b 2 , the real emission term with a factor 
Jo(|fcj_||6|) cancels the virtual emission term with a factor 1 because Jo(|fej_||6|) — 1 ~ 
0, as in our other analyses. Thus this region does not contribute to the evolution of 
Q(b, Y; ?7a, rjb, a, b). In the region k\ S> 1/b 2 , only virtual emissions contribute to the 
evolution of Q(b, Y; ry a , r/b, a, b). The rapidly oscillating Jo(|fcj_||6|) factor that multiplies 
Tlf cip eliminates the real emission contribution. In the region k"j_ ~ 1/b 2 , both real and 
virtual emissions contribute. 

We have by no means carried out a detailed investigation of this angle ordered scheme. 
However, it appears to us that it will give equivalent results for the Z-boson transverse 
momentum distribution as does the virtuality ordered shower. 
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13.6 Dipole antenna shower 

In this paper so far we have discussed dipole parton showers in which the creation of a 
new gluon is attributed to the splitting of one of the previously existing partons. There is 
an ambiguity because, for soft gluon emissions, interference diagrams are important: one 
has to include emission from a parton I squared, emission from a parton k squared, and 
the l-k interference graphs. Using the eikonal approximation, the sum of these is simple, 
as in eq. (13.1). We have partitioned the total emission probability into a fraction A' lk 
associated with splitting of parton I and a fraction A' kl = 1 — A' lk associated with the 
splitting of parton k, as in eq. (13.2). In the nomenclature of section 16 of ref. [34], this 
constitutes a partitioned dipole shower. 18 There is a separate momentum mapping Vi for 
each parton I that splits. Here V\ is an operator, defined in ref. [12], on the space of parton 
states. It is not necessary that this momentum mapping depends on the choice of k, and 
the momentum mapping used in this paper and ref. [12] does not depend on k. Thus, 
in a highly compressed notation, we can write that the emission from the l-k dipole is 
partitioned into two terms, 

Hf* rt (t) oc Vi A\ k + V k A' kl ^ . (13.19) 

Pm+l 'Pi Pm+1 -Pk Pm+1 'Pi Pm+1 'Pk 

In section 13.2, we have seen that the summation of large logarithms can be sensitive 
to the choice of the partitioning function A' lk ; a choice based on the angles between parton 
I and k and the emitted soft gluon is preferred. 

One might be tempted to get rid of this ambiguity, particularly in the leading color 
approximation where one can use pairs of color connected partons as the dipoles. The idea 
is to consider each l-k dipole as a unit that can emit a gluon. Thus the basic building 
blocks are 2 — ► 3 parton splittings. A shower based on this approach may be called a dipole 
antenna shower. The pioneering development along these lines is the final state shower of 
Ariadne [35]. More recent examples include those in ref. [19]. There is a corresponding 
subtraction scheme for next-to-leading order calculations, antenna subtraction [36]. 

In a dipole antenna shower, there is no A' lk . There is a separate momentum mapping 
Vik for each dipole l-k. Thus eq. (13.19) is replaced by 

Hff(t) oc V ik * ^ • (13.20) 

Pm+1' Pi Pm+l-Pk 

The freedom to choose A' lk now resides in the freedom to choose Vik- It must be symmetric 
under the interchange of partons / and k. It is usually defined in such a way that momentum 
is conserved locally in the 2^3 splitting, without taking momentum from any of the other 
partons. For a final state dipole, this means 

Pl+Pk=Pl+Pk+Pm+l • (13.21) 

Such a mapping is defined in ref. [37]. We note that it is possible to have a shower that 
is simultaneously a partitioned dipole shower and an antenna dipole shower. We get that 



18 Note that from the point of view adopted in section 13.5, an angle ordered shower also can be considered 
as a partitioned dipole shower. 
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case if we define 

Viu = 0(<V+i < ^ + 0($ k , m+1 < <?, jm+1 ) P fc , (13.22) 

where the is the angle between momenta i and j and 7^ represents the momentum 
mapping that is used in this paper. Then A' lk = 0(i9z )m +i < ~&k,m+i)- 

As far as we can see, a dipole antenna shower can reproduce the proper summation of 
large logarithms for the Z-boson transverse momentum distribution studied in this paper. 
However, the momentum mapping for initial state splittings cannot be local in the sense 
of eq. (13.21). For a dipole consisting of the two initial state partons, eq. (13.21) becomes 

Pa +Pb = Pa + Pb ~ Pm+l ■ (13.23) 

We want p^± = p^± = and also j?b,± = Pb,± = 0. This leaves no place for the transverse 
momentum of parton m + 1 to go. One might hope to use a non-local momentum mapping 
for this case, but use a local mapping in the case of a dipole consisting of one initial state 
parton, say "a", and a final state colored parton k. In this case, eq. (13.21) becomes 

Pa~ Pk = Pa~ Pk ~ Pm+l ■ (13.24) 

However, that would mean that pk,± — Pk,± = ~Pm+i,i : the recoil transverse momentum 
from the emission of parton m + 1 is taken up by parton k instead of by the Z-boson. 
Thus a non-local momentum mapping is needed for all dipoles that include an initial state 
parton. 



14. Conclusions 

In this paper, we have examined the differential cross section da / {dp+dY) for producing 
a Z-boson in hadron-hadron collisions as a function of the transverse momentum p± and 
rapidity Y of the Z-boson. For p\ <^ M|, the perturbative expansion of da/(dp±dY) 
contains large logarithms, log(p5_/M|). There are known QCD results for the summation 
of these logarithms, quoted in section 12. We have asked to what extent a virtuality ordered 
parton shower algorithm of the sort defined in ref. [12] (with some small modifications as 
discussed in this paper) correctly reproduces the known results. 

The parton shower evolves in shower time t, defined to be minus the logarithm of the 
virtuality of a parton splitting, so that hard splittings come first, soft splittings last. The 
state of the shower at time t, in the sense of the distribution of parton configurations, is 
represented by a vector |/o(t)). In this notation, da / (dp±dY) as it has developed by shower 
time t is denoted by (p±, y|p(t)). Then, at the end of the shower at time if, the predicted 
cross section is 

We have examined whether the cross section obtained from the shower algorithm has the 
structure of the QCD result by examining how (p_i_, y|/o(i)) evolves with t. 

To find d(pj_, Y \p(t)) /dt, we used eq. (2.9), which expresses how the completely exclu- 
sive partonic state evolves. This expresses d(p±, Y\p(t)^j /dt in terms of the complete state 
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of all of the partons at time t. Fortunately, we found that, with appropriate approxima- 
tions, d(p±, Y \p(tf)/dt is expressed as a convolution of a certain kernel with (p±, Y . 
That is, the detailed configuration of all of the partons does not matter, only the inclu- 
sive distribution (p±,Y\p(t)"j matters. Thus we obtain a closed form differential-integral 
equation for [p±, Y\p(t)^j . 

The result can be simply stated in terms of the Fourier transform of transverse mo- 
mentum distribution, defined in eq. (4.4), 

(l\Q(b,Y)\p(t)) e- ib - p -(p±,Y\p(t)) . (14.2) 

We found that (l | Q(b, Y) p(i)) stops evolving at a certain evolution time t c , with t c = 
log (6 2 M 2 e 27E /4). We found in eq. (12.1) that (l\Q(b,Y)\p(t c )) is a convolution of parton 
distribution functions evaluated at the appropriate scale times a function that does not 
involve parton distribution functions, 



1 Q(b,Y )\p (t c )) = / *fc / *fc U/AU ~ A , , J, — 

~^ Jo Jo 4n c (a)n c (b)2ri a r] h p A -pB 



(14.3) 



(l\Q(b,Y;r] a ,r} h ,a, b)\p pcrt (t c )) . 



The important structure is thus expressed in the function (l| Q(b, Y; rj a , rjb, a, b)\p peT t(t c )) . 

Consider the logarithm of (l| Q(b, Y; rf a , rjb, a, b) Lo per t(ic)) • On general grounds, one 
would expect that this function has two powers of the large logarithm log(Af 2 b 2 ) per power 
of a s (M 2 ): 



log[(l| Q(b, Y; ?7 a ,r? b , a, b)\p pelt (t c ))] 



m— 1 m— N / 



n=l m=0 

However, we found that the partonic cross section produced by the shower algorithm ex- 
ponentiates in 6-space in the sense that 

D nm = for m > n + 1 . (14.5) 

This property is shared with the full QCD result. 

This exponentiation is not guaranteed by simply having a parton splitting probability 
that matches QCD in the soft and collinear limits, including the proper interference between 
emissions of a soft gluon from different partons. We found that two other features of the 
shower algorithm are important. 

First, partons in a parton shower are treated as being on-shell, but it is not possible for 
an on-shell parton to split into two on-shell partons and still conserve momentum. There 
has to be a momentum mapping that takes a small amount of momentum from somewhere 
and supplies it to the daughter partons. In the case of an initial state splitting, the problem 
is more pronounced because the initial state partons are treated as having zero transverse 
momentum. This means that when a soft gluon is emitted from an initial state parton, 
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some transverse momentum has to come from somewhere. We found that the momentum 
mapping has to be such that the Z-boson gets most of the recoil transverse momentum. 
In fact, the momentum mapping proposed in ref. [12] did not have this property and 
thus needed to be modified. Similarly, we noted that the momentum mapping used for 
perturbative subtractions in the Catani-Seymour dipole subtraction scheme does not give 
the desired exponentiation for (l|Q(b, Y)\p(t c )^ in a dipole based shower. 

Second, in a shower like that of ref. [12] that is based on color dipoles, there is a 
function that we call A' lk that specifies how much of the radiation from dipole l-k is treated 
as being emitted from parton I and how much is treated as being emitted from parton k. 
We found that, for initial state emissions, the choice of A' lk matters. A choice based on the 
angles between the momenta of partons I and k and an emitted soft gluon is satisfactory, 
while the choice used in the Catani-Seymour dipole subtraction scheme fails to give the 
proper exponentiation in a dipole based shower. 

After having checked for the proper exponentiation of the large logarithms, the next 
question was which non-zero coefficients D nm are correctly produced by the shower algo- 
rithm. 

We found that the leading coefficient, -D12, is correctly reproduced. This result has a 
simple physical interpretation in the parton shower picture. The quantity (l | Q(b, Y) \p{t{ )) 
is the cross section to produce the Z-boson with rapidity Y but not radiate initial state 
gluons with transverse momentum bigger than approximately 1 /b 2 . The Sudakov exponent 
S is the integral of the differential probability to radiate such a gluon, so exp(— S) is the 
probability not to radiate any such gluons. The leading term in the Sudakov exponent is 
Cpa s /TT times twice the area of the triangle in figure 2. This area is total phase space for 
gluon emission with k 2 ^ > 1/b 2 from one of the two initial state partons; doubling the area 
gives the total phase space for emission from either line. Quantum coherence plays a role 
here. The initial emission is from a color dipole consisting of the two incoming quarks. 
Emissions with large positive rapidity count as emission from incoming quark "a" while 
emissions with large negative rapidity count as emission from incoming quark "b" . 

We found also that the next-to-leading log coefficient, D±\, is correctly reproduced. 
This coefficient corresponds to the behavior of the parton splitting near the three edges 
of the triangle in figure 2. It thus represents physics that is more subtle than the physics 
behind the leading coefficient, D\2- 

We note that once the choices of momentum mapping and A' lk have been made, no 
adjustment of parameters is needed to get D\i and Dn to match the QCD result. 

One could argue that matching the leading order coefficients Dyi and D\\ is all that 
one should expect from a leading order parton shower. Nevertheless, we found that it 
is possible to do better (as in ref. [7]). With a suitable choice of the argument of a s in 
the splitting function, it is possible to match all of the coefficients D nm with m = n + 1 
and with m = n. With this level of matching, the Sudakov exponent S produced by the 
shower algorithm is a good approximation to the full QCD exponent Sqcd m t ne limit 
log(M 2 6 2 ) 00 and a s (M 2 ) -> with a s (M 2 ) log(M 2 fo 2 ) fixed. 

We have investigated how well the virtuality ordered parton shower algorithm defined 
in refs. [12, 13, 14] and this paper performs in summing large logarithms for the transverse 
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momentum of a vector boson produced in the Drell-Yan process in hadron-hadron colli- 
sions. The same analysis would apply to the transverse momentum distribution of Higgs 
bosons produced in hadron-hadron collisions. There are other sorts of summations of large 
logarithms that are of importance for understanding experiments. We believe that the in- 
vestigation of other cases, perhaps using methods developed in this paper, is an important 
goal, both for the virtuality ordered parton shower algorithm used in this paper and for 
other parton shower algorithms. 

Acknowledgments 

This work was supported in part by the United States Department of Energy and the 
Helmoltz Alliance "Physics at the Terascale." 

A. Structure of inclusive splitting 

In this appendix, we review formulas from refs. [12, 13, 14] for inclusive splitting, that is 
the splitting operator l H\{t) times the inclusive measurement function (l|. We start with 
in eq. (12.20) of ref. [12], 

(l\Hi(t)\{p,f,i/,<f,8,c} m ) =2({s'} m \{s} m )({c'} m \h(t,{p,f} m )\{c} m ) . (A.l) 

The inclusive splitting function is diagonal in spin, but it has a non-trivial color structure. 
The structure of the splitting is contained in the operator h, which is given in eq. (12.21) 
of ref. [12]. With a slightly modified notation and with one error corrected, this equation 
reads 

({c'} m \h(t,{p,f} m }\{c} m } = 



2E E / d C P #(C P e/XM m ,Cf)) 

Ql 



2 

l Cf6*i(/i) 



x5\t-\o 



x 



|(p, + (-l)*«..+*«.'.p ro+1 )2- m 2(/ l )| 

{a)n c (b) 77 a 77b vVih/B^ Vf) 



(A.2) 



n c (a)n c (b) fa/A(v^ ^h/B^b, Mf) 

X j#(/m+l /g) ({c'}m\{c} m ) C(fl, f m+ i) W U ({f,p} m+1 ) 

- 6(f m+1 =g)J2 ({c'}m\T k • 7]|{ C } m ) 

k^l 

X [wil{{f,p}rn+l) ~ 2Ai k ({p} m+1 )wik({f,p} m+ l)] 

Here possible quark masses are included. There is a sum over the index I of the initial state 
or final state parton that splits. There is a sum over flavor choices £f for the splitting and an 
integration over the momentum choices, J d(, p . Once we choose splitting variables t, z and 
(j), this becomes an integration over the splitting variables. There is a delta function that 
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expresses the definition of the shower time t. Next, there is a ratio of parton luminosities, 
containing parton distribution functions, as in eq. (4.8). The main momentum and color 
dependence is contained inside the braces that follow. 

The second term inside the braces applies when the newly emitted parton with label 
m + 1 is a gluon. There is a sum over partons k, where k ^ I. Partons I and k constitute 
a color dipole that can emit the gluon. There are functions wu and wik that describe the 
direct terms and the interference terms, respectively, for gluon emission. These functions 
are defined in ref. [12]. Convenient expressions for wu and wik as rational functions of dot 
products of the momenta involved are given in ref. [13]. 19 There is also a function Aik 
that expresses how the interference term is partitioned between a part considered to be 
a splitting of parton I and a corresponding part considered to be a splitting of parton k. 
There is a color operator T k • 7) that inserts a gluon color matrix T a on line k and another 
on line I and sums over the gluon color index a. This operator was denoted gik({f}m+l) 
in the original equation. 

The first term inside the braces applies when the newly emitted parton with label 
m + 1 is not gluon. Then there is no interference term. There is a color factor 



C(fl, fm+l) 



Tr {fi, /m+i} = {<?, q} or {q, q} 

C F {f l J m+1 } = {g,q},{g,q},{q,g}or{q,g} . (A.3) 

Ca {fl,f m+ l} = {g,g} 



The cases with f rn +i = g do not appear in eq. (A. 2), but we will need these cases later. 
In eq. (12.21) of ref. [12], we failed to note the possibility of an initial state splitting with 
{//, fm+i] = {g, q] or {g, q} } so we listed this factor simply as T R . 

We now manipulate eq. (A. 2) so as to obtain a more useful form. For the f m +i = 
g term, we divide the splitting function into two parts by adding and subtracting the 
eikonal approximation to wu. (The eikonal approximation is the limiting form when the 
gluon momentum tends to zero. The function wik is already constructed in the eikonal 
approximation in ref. [12].) Thus we write 

wii{{f,p} m+ i)-2Ai k ({p} m+ i)wi k ({f,p} m+ i) 

= (ia({/,pWi) - ^/ konal ({/,pWi)) (A.4) 
+ (™f onal ({/,pWi) - 2A lk ({p} m+1 )w lk ({f,p} m+1 )) . 

The first term is important in only in the collinear limit. The momentum dependent 
factor w tt - wjf konal is independent of k and multiplies ({c'} m |Tfc • Tj|{c} m ). When we sum 
over k, we can use 

{{c}m\T k ■ D|{c} m } = - ({c'} m \Ti • 7]|{c} m ) 
W (A.5) 

= -C{{f u f m+l }) ({c'} m \{c} m ) . 



9 In ref. [13] and also in ref. [14], we use the notations Wu = wu and Wik = 2^4; 



\kWlk- 
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For the second term, we take Aij. to be of the form specified in eq. (7.2) of ref. [14], in 
which Aik is expressed using another function A',, . This gives the result of in eq. (7.10) of 
ref. [14], 

^ konal ({/,p}m+i)-2^ fc ({p} m+1 )^ fc ({/,p} m+1 ) 

-PI ( A - 6 ) 



4vra s A' lk ({p} m+1 ) 



(Pm+1-Pl Pm+l-Pk 



1 2 



where 



D lk = Pm+l 'PI Pk ~ Pm+l -Pk Pi ■ (A. 7) 



This is for general masses. For the analysis of the Z-boson p± distribution, we want all of 
the partons to be massless. Then 

(^f^({/,p} m+1 )-2A^({p} m+1 )^ fc ({/,p} m+1 )) 

= 4-!ra s A lk {{p} m+1 ) — — . 

p m+1 -pi p m+ i -p k 

Our choice in this paper for the function A' lk is given in eq. (13.4). 
With these rearrangements, we have 



(l\Hi{t)\{p, /, s,c,s,c} rn ) 

= ({s'} m \{s} m ) EE /^P G r KMm,Cf)) 



Cf6#i(/i) 

x S [ t - lo 



Ql 



|(K + (-l)^+^p m+1 )2_ m 2 (/i) | 

n c (a)n c (b) 1]eJ]h fa/A^a, ^ i2 F)fb/B^ h ^• ^f) 
n c (a)n c (b) fj^ fa/A^m M|)/b/s( 7 ?b' /40 

X j({ c '}m|{c}m} C(fi,f m +l) 

x W«({/,p}m+l) - 0(L+l = g)wff°™ l ({f,p} m+1 ) 
- 6(f m+1 = g) E <{c'} m |T fc • Ti\{c} m ) 

k^l 

-P 2 

x 47ra B A{ fc ({p} m+ i) 



(A.9) 



{Pm+l 'Pi Pm+l 'Pk 



\ 2 



We can now specialize to initial state splittings with massless partons. In the sum over 
the label I of the parton that splits, we take I = a. The integration measure d^ p , with the 
choice of variables {t, z, <p} used in this paper, is 

d( p = 2p a • p b (27r) dtdz — — = — — dtdz — . (A.1U) 

27T 4z 07T Z 2 2"7T 

The integration over t is performed using the delta function that defines t. We have de- 
fined the matrix element of 7ia CTt (t, z, 4>, /') in eqs. (4.15) and (6.2) to be the coefficient of 
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dz d(p / (2tt) in the integrand and sum over flavors, leaving out the ratio of parton luminosi- 
ties. Cf. eqs. (8.4) and (8.5). Thus 



Pm+1 'Pa 
8vr 2 



^{s }m | {^}m} 



X |({c / }m|{c} m )C(/ a ,/ m+ i) 

X W aSL {{f,p} m+ l) ~ 9{f m+ i = g) 

-9(f' = g)Y,({c} m \T k -T a \{c} m ) 

fc^a 

1p~aPk 



cikonal 



({f,p} 



m+1, 



(A.ll) 



x 4vra s A' ak ({p} m+1 ) 

p m+ i -p a p m+1 -p k 

We can use the results of ref. [13] to state the functions «J aa in terms of the variables 
y,z,(f> used in this paper. For {/ a ,/m+i} = {<7,g} or {q,g}, we find, using eq. (2.29) of 
ref. [13], 



Pm+1- Pa, 

8tt 2 



:ikonal 



aa "'aa 



Os 1 

2vr z 



(l-z)(l + y) 



For {/a,/m.+i} = {g,g}, we find, using eq. (2.57) of ref. [13], 

as 2 
2vr z 



Pm+VPa 

8vr 2 



w aa - ™ aa 



cikonal 



— + <i--)- 7 f fi ^- 

2 [l-z) + y 



For {/ a ,/m.+i} = {<?, q} or {g, g}, we find, using eq. (A.2) of ref. [13], 

q s 1 l + (l-z) 2 



Pm+1 'Pa, — 

8vT 2 W °* ~ 27T Z 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



For {fa,fm+i} = {g,q} or {g,q}, we find, using eq. (A.3) of ref. [13], 

^^aa = ^ ^ [, 2 + (1 - ,) 2 ] • 

8ir z 2-ir z 

The contribution to (l|zW a ert (t, z, (ft, f')\{p, /, s' , c' , s, c} m ) from the term proportional 
to A' ak in the last line of eq. (A.ll) vanishes if the emitted parton is not a gluon and if 
/' = g it is 



[l\znr\t,z, 0,g)|{p, f,s', d, s, c] m ) 



cikonal 



{{s'}m\{s}m)(-l) £ <{C%|T, • T a |{ C } m > £ < fc ({p} m+1 ) 



(A.16) 



fc^a 



Pm+1 -Pfc 



Let r k be the rapidity of parton k in the p a + p^ rest frame and let (pk its azimuthal angle, 
while 4> denotes the azimuthal angle of parton m + 1. Then, using the definition (13.4) of 
A' lh from eq. (7.12) of ref. [14], we find 

{l\z-Hr\t, z, </>, g)\{p, f, s', c', s, c} m ) eikonal 

= ({s'}m\{s}m)(-l) ({ C '}m\Tk ■ T a \{ C } r 

fc^a 



2-7T 1 — z + y 



f{z,y,(p,r k ) 



(A.17) 
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where 

1 _ 2y/(l- z)y 



f(z,y,<f>,r k ) 



2 

cos(0 - (j) k ) 



1 + y 1 — z + y 1 + y 

i ( A -!8) 



2y z 2 



(1-z + y) (1 + y) 2 



In the limit that the emitted gluon is soft, we can neglect (1 — z) compared to 1 and 
y compared to 1. We do not, however, neglect y compared to (1 — z). Then 



f(z,y,(/>,r k ) 



1 - e Tk — f- ■ cos(0 - 4> k ) + e 



1-z + y ' ' (1-z + y) 



-i 



(A.19) 



We use this result in section 7.1. 

When the emitted gluon becomes collinear, y — > with fixed z, we get 

f(z,0,<P,r k ) = l . (A.20) 

Then we can use eq. (A. 5), so that 

{l\znr\t, z, 0, g)|{p, /, a', c', s, c} m ) eikonal 
~ (W}m\{s}m){{c'}m\{c} m ) C(fi, f m +l) ^ - x _^ + y ■ 

The behavior of the splitting probability in the collinear limit is obtained by combining 
this result with the results in eqs. (A. 12), (A.13), (A.14), and (A. 15). 

B. An integral of Bessel functions 

In this appendix, we prove a theorem about the Jo Bessel function that is needed in 
section 9.3. The theorem concerns an integral of the form 



F(x 2 , Xl ) = r — /(logs) [J (x) - 1] . (B.l) 



where x here stands for |fcj_||b|. We imagine that ii < 1 and x 2 ^> 1. We suppose 
that the function /(log a;) is slowly varying in the sense that its second derivative is small: 
/"(log a;) <C /(log a;). In our application, / is o s (fc 2 L ) = a s (x 2 /b 2 ). This is slowly varying 
because its second derivative with respect to log(x) is proportional to a^. If the integration 
extends down to small values of k 2 ^, then a s (fc 2 L ) is not slowly varying if we use the 
perturbative evolution equation, so the theorem is useful only if we use an approximation 
in which a s (fc 2 L ) stops varying for small values of k 2 ^_. 

The theorem concerns the approximation of F(x 2 , xi) by 

f X2 dx 

i ? (x2,xi) a pp rox = - / — /(log 9(x > xq) , (B.2) 

Jxi x 

where 

x = 2e _7E . (B.3) 
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We thus analyze the difference, 

AF(x 2 ,x\) = F(x 2 ,xi) - F(x2,xi) appiox , (B.4) 

which is given by 

AF(x 2 ,xi)= — f(\ogx)[J o (x)-0(x<x o )} . (B.5) 

J xi x 

We expect that AF(x2, x\) is small because for the integration range Jo(x) oscillates, 

so that its integral is small, while for the integration range iCl, Jo(%) ~ 1> so that it 
cancels the term 6(x < xq). 

To prove that AF(x2,x\) is small, we define 

r x df 

h{x)= / — [J (x)-9(x<x )} . (B.6) 
Jo x 

Now, I\(x) vanishes as x — > 0. The approach to the limit is given by 

x 2 

h(x) ~ -— , x -> . (B.7) 



The integral converges for x — * oo and, with our special choice of xq, its value is Ji(oo) = 0. 
The approach to the limit is given by 

h(x) ~ -\ ^sin - x) , x -> oo . (B.8) 



7T £ 3 / 2 \4 

We rewrite AF(x2, x±) using an integration by parts: 



f X2 dx d 
AF(x 2 ,xi)= J — /(logx)x— h(x) 



[f(logX2)h(x2)-f(logx 1 )h(x 1 )} (B.9) 

x ' 2 dx 

— /'(log a) Ji(x) . 

Now define 

h(x) = / — log(x/x) [J (x) - 0(x < x )] . (B.10) 
Jo x 

Note that the derivative of I 2 (x) with respect to logx is I\(x). Again, /2(0) = 0. The 
approach to the limit is given by 

9 

I 2 (x)~~ , x->0 . (B.11) 
lb 

Again, the integral converges for x — > oo and its value is, perhaps surprisingly, 12(00) = 0. 
This can be proved using the generating function in eq. (4.12) of ref. [38]. The approach 
to the limit is given by 

h(x) ~ -\ fts-cos - x) , x ^ 00 . (B.12) 



7T X 5 / 2 V4 
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We can write AF{x 2 , x\) using 



AF(x 2 ,x 1 ) 



[/(log x 2 ) h (x 2 ) - /(logxi)Ji(xi)] 



r x 2 dx , d 
-/ — f'(logx)x—I 2 (x) 
J X1 x dx 



(B.13) 



This gives 



AF(x 2 ,x 1 ) 



[f(logx 2 )h(x 2 ) - /(logxi)Ji(xi)] 
- [/'(logx 2 ) h{x 2 ) - f(\ogx\) h{xi)] 

+ / - f"{\ogx)I 2 {x) . 



(B.14) 



Now we can see the conditions under which AF(x 2 ,xi) is small. There are terms 
proportional to /(loga^) and /'(log x 2 ). Using the large x behavior of I\(x 2 ) and I 2 (x 2 ), 
we see that these contributions are power suppressed for large x 2 as long as /(logX2) and 
/'(log x 2 ) do not grow for large x 2 . There are terms proportional to /(log x\) and /'(log x\) . 
Using the small x behavior of Ii{x\) and -^2 (^1 ) , we see that these contributions are power 
suppressed for small x\ as long as /(logxi) and /'(log^i) do not grow for small x\. Finally, 
there is an integral of /"(log x) I 2 {x). Since I 2 (x) falls off for large x and for small x, it 
is the behavior of /"(logx) around x = 1 that is important in the integration (as long as 
/"(logx) is bounded everywhere). If /"(logx) is small, AF(x 2 , xi) will be small. 

Additionally, AF(x 2 ,xi) is small when both x\ and x 2 are large compared to 1. In 
this case, all of the terms in eq. (B.14) are suppressed by powers of 1/xi and l/x2- 
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